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Recently, it was pointed out that quantum orders and the associated projective symmetry groups 
can produce and protect massless gauge bosons and massless fermions in local bosonic models. In 
this paper, we demonstrate that a state with such kind of quantum orders can be viewed as a 
string-net condensed state. The emerging gauge bosons and fermions in local bosonic models can 
be regarded as a direct consequence of string-net condensation. The gauge bosons are fluctuations 
of large closed string-nets which are condensed in the ground state. The ends of open strings (or 
nodes of open string-nets) are the charged particles of the corresponding gauge field. For certain 
types of strings, the nodes of string-nets can even be fermions. According to the string-net picture, 
fermions always carry gauge charges. This suggests the existence of a new discrete gauge field that 
couples to neutrinos and neutrons. We also discuss how chiral symmetry that protects massless 
Dirac fermions can emerge from the projective symmetry of quantum order. 

PACS numbers: 11.15.-q 



Contents 



I. Introduction 1 

A. Fundamental questions about light and 
fermions 1 

B. Gapless phonon and symmetry breaking 
orders 2 

C. The existence of light and fermions impUes 
the existence of new orders 2 

D. Topological order and quantum order 2 

E. The quantum orders from string-net 
condensations 3 

F. Organization 5 

II. Local bosonic models 5 

III. Z2 spin liquid and string-net 
condensation on square lattice 6 

A. Hamiltonians with closed-string-net 
condensation 6 

B. String condensation and low energy effective 
theory 7 

C. Three types of strings and emerging 
fermions 7 

IV. Classification of diff"erent string 
condensations by PSG 8 

A. Four classes of string-net condensations 8 

B. PSG and ends of condensed strings 9 

C. PSG's classify different string-net 
condensations 10 

D. Different PSG's from the ends of different 
condensed strings 11 

V. Massless fermion and PSG in string-net 
condensed state 11 

A. Exact soluble spin-^i model 12 



*URL: http: //dao .mit . edu/'wen 



B. Physical properties of the spin-ii model 13 

C. Projective symmetry and massless fermions 14 

VI. Massless fermions and string-net 

condensation on cubic lattice 16 

VII. Artificial light and artificial massless 

electron on cubic lattice 16 

A. 3D rotor model and artificial light 16 

B. (Quasi-)cxact soluble QED on cubic lattice 17 

C. Emerging chiral symmetry from PSG 18 

VIII. QED and QCD from a bosonic model on 

cubic lattice 20 

IX. Conclusion 22 

References 23 



I. INTRODUCTION 

A. Fundamental questions about light and fermions 

We have known light and fermions for many years. But 
we still cannot give a satisfactory answer to the following 
fundamental questions: What are light and fermions? 
Where light and fermions come from? Why light and 
fermions exist? At moment, the standard answers to the 
above fundamental questions appear to be "light is the 
particle described by a gauge field" and "fermions are 
the particles described by anti-commuting fields" . Here, 
we like to argue that there is another possible answer to 
the above questions: our vacuum is filled with string- like 
objects that form network of arbitrary sizes and those 
string-nets form a quantum condensed state. Accord- 
ing to the string-net picture, the light (and other gauge 
bosons) is a vibration of the condensed string-nets and 
fermions are ends of strings (or nodes of string- nets). The 
string-net condensation provides a unified origin of light 
and fermions. [62] 
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Before discussing the above fundamental questions in 
more detail, we would like to clarify what do we mean by 
"light exists" and "fermions exist" . We know that there 
is a natural mass scale in physics - the Planck mass. 
Planck mass is so large that any observed particle have a 
mass at least factor 10^^ smaller than the Planck mass. 
So all the observed particles can be treated as massless 
when compared with Planck mass. When we ask why 
some particles exist, we really ask why those particles are 
massless (or nearly massless when compared with Planck 
mass). So the real issue is to understand what makes 
certain excitations (such as light and fermions) massless. 
We have known that symmetry breaking is a way to get 
gapless bosonic excitations. We will see that string-net 
condensation is another way to get gapless excitations. 
However, string- net condensations can generate massless 
gauge bosons and massless fermions. 

Second, we would like to clarify what do we mean by 
"origin of light and fermions" . We know that everything 
has to come from something. So when we ask "where 
light and fermions come from", we have assumed that 
there are some things simpler and more fundamental than 
light and fermions. In the section II. we define local 
bosonic models which are simpler than models with gauge 
fields coupled to fermions. We will regard local bosonic 
models as more fundamental (the locality principle). We 
will show that light and fermions can emerge from a local 
bosonic model if the model contains a condensation of 
nets of string-like object in its ground state. 

After the above two clarifications, we can state more 
precisely the meaning of "string-net condensation pro- 
vides another possible answer to the fundamental ques- 
tions about light and fermions". When we say gauge 
bosons and fermions originate from string-net condensa- 
tion, we really mean that (nearly) massless gauge bosons 
and fermions originate from string-net condensation in a 
local bosonic model. 



B. Gapless phonon and symmetry breaking orders 

Before considering the origin of massless photon and 
massless fermions, let us consider a simpler massless (or 
gapless) excitation - phonon. We can ask three simi- 
lar questions about phonon: What is phonon? Where 
phonon comes from? Why phonon exists? We know that 
those are scientific questions and we know their answers. 
Phonon is a vibration of a crystal. Phonon comes from 
a spontaneous translation symmetry breaking. Phonon 
exists because the translation-symmetry-breaking phase 
actually exists in nature. In particular, the gaplessness of 
phonon is directly originated from and protected by the 
spontaneous translation symmetry breaking. [1, 2] Many 
other gapless excitations, such as spin wave, superfluid 
mode etc , also come from condensation of point-like ob- 
jects that break certain symmetries. 

It is quite interesting to see that our understanding of a 
gapless excitation - phonon - is rooted in our understand- 
ing of phases of matter. According to Landau's theory, [3] 
phases of matter are different because they have different 
broken symmetries. The symmetry description of phases 



is very powerful. It allows us to classify all possible crys- 
tals. It also provides the origin for gapless phonons and 
many other gapless excitations. Until a few years ago, it 
was believed that the condensations of point-like objects, 
and the related symmetry breaking and order parame- 
ters, can describe all the orders (or phases) in nature. 

C. The existence of light and fermions implies the 
existence of new orders 

Knowing light as a massless excitation, one may won- 
der maybe light, just like phonon, is also a Nambu- 
Goldstone mode from a broken symmetry. However, ex- 
periments tell us that a U{1) gauge boson, such as light, 
is really different from a Nambu-Goldstone mode in 3-1-1 
dimensions. Therefore it is impossible to use Landau's 
symmetry breaking theory and condensation of point- 
like objects to understand the origin and the massless- 
ness of light. Also, Nambu-Goldstonc modes are always 
bosonic, thus it is impossible to use symmetry breaking 
to understand the origin and the (nearly) masslessness 
of fermions. It seems that there does not exist any order 
that can give rise to massless light and massless fermions. 
Because of this, we put light and electron into a different 
category than phonon. We regarded them as elementary 
and introduced them by hand into our theory of nature. 

However, if we believe light and electrons, just like 
phonon, exist for a reason, then such a reason must be a 
certain order in our vacuum that protect the masslessness 
of light and electron. (Here we have assumed that light 
and electron are not something that we place in an empty 
vacuum. Our vacuum is more like an "ocean" which is 
not empty. Light and electron are collective excitations 
that correspond to certain patterns of "water" motion.) 
Now the question is that what kind of order can give rise 
to light and electron, and protect their masslessness. 

If we really believe in the equality between light, elec- 
tron and phonon, then the very existence of light and 
fermions indicates that our understanding of states of 
matter is incomplete. We should deepen and expand our 
understanding of the states of matter. There should be 
new states of matter that contain new kind of orders. The 
new orders will produce light and electron, and protect 
their masslessness. 



D. Topological order and quantum order 

After the discovery of fractional quantum Hall (FQH) 
effect, [4, 5] it became clear that the Landau's symmetry 
breaking theory cannot describe different FQH states, 
since those; states all have the same symmetry. It was 
proposed that FQH states contain a new kind of order - 
topological order. [6] Topological order is new because it 
cannot be described by symmetry breaking, long range 
correlation, and local order parameters. Non of the usual 
tools that we used to characterize phases applies to topo- 
logical order. Despite of this, topological order is not an 
empty concept. Topological order can be characterized 
by a new set of tools, such as the number of degenerate 
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FIG. 1: A classification of different orders in matter. (We 
view our vacuum as one kind of matter.) 



ground states, qiiasiparticlc statistics, and edge states. It 
was shown that the ground state degeneracy of a topo- 
logical ordered state is a universal property since the 
degeneracy is robust against any perturbations. [7] Such 
a topological degeneracy demonstrates the existence of 
topological order. It can also be used to perform fault 
tolerant quantum computations. [8] 

Recently, the concept of topological order was general- 
ized to quantum order. [9, 10] Quantum order is used to 
describe new kinds of orders in gapless quantum states. 
One way to understand quantum order is to see how it 
fits into a general classification scheme of orders (sec Fig. 
1). First, different orders can be divided into two classes: 
symmetry breaking orders and non-symmetry breaking 
orders. The symmetry breaking orders can be described 
by a local order parameter and can be said to contain 
a condensation of point-like objects. All the symmetry 
breaking orders can be understood in terms of Landau's 
symmetry breaking theory. The non-symmetry breaking 
orders cannot be described by symmetry breaking, nei- 
ther by the related local order parameters and long range 
correlations. Thus they are a new kind of orders. If a 
quantum system (a state at zero temperature) contains 
a non-symmetry breaking order, then the system is said 
to contain a non-trivial quantum order. We see that a 
quantum order is simply a non-symmetry breaking order 
in a quantum system. 

Quantum order can be further divided into many sub- 
classes. If a quantum state is gapped, then the corre- 
sponding quantum order will be called topological or- 
der. The low energy effective theory of a topological 
ordered state will be a topological field theory. [11] The 
second class of quantum orders appear in Fermi liquids 
(or free fermion systems). The different quantum or- 
ders in Fermi liquids are classified by the Fermi surface 
topology. [10, 12] 



E. The quantum orders from string-net 
condensations 



In this paper, we will concentrate on the third class of 
quantum orders - the quantum orders from condensation 
of nets of strings, or simply, string- net condensation. [13, 
14] This class of quantum orders shares some similar- 
ities with the symmetry breaking orders of "particle" 
condensation. We know that different symmetry break- 



ing orders can be classified by symmetry groups. Using 
group theory, we can classify all the 230 crystal orders 
in three dimensions. The symmetry also produces and 
protects gapless Nambu-Goldstone bosons. Similarly, as 
we will see later in this paper, different string-net con- 
densations (and the corresponding quantum orders) can 
be classified by a mathematical object called projective 
symmetry group (PSG).[9, 10] Using PSG, we can clas- 
sify over 100 different 2D spin liquids that all have the 
same symmetry. [9] Just like symmetry group, PSG can 
also produce and protect gapless excitations. However, 
unlike symmetry group, PSG produces and protects gap- 
less gauge bosons and gapless fermions.[9, 15, 16] Because 
of this, we can say light and massless fermions can have 
a unified origin. They can come from string-net conden- 
sations. 

We used to believe that to have light and fermions in 
our theory, we have to introduce by hand a fundamen- 
tal J7(l) gauge field and anti-commuting fermion fields, 
since at that time we did not know any collective modes 
that behave like gauge bosons and fermions. Now, we 
know that gauge bosons and fermions appear commonly 
and naturally in quantum ordered states, as fluctuations 
of condensed string-nets and ends of open strings. This 
raises an issue: do light and fermions come from a fun- 
damental ?7(1) gauge field and anti-commuting fields as 
in the 123 standard model or do they come from a par- 
ticular quantum order in our vacuum? Clearly it is more 
natural to assume light and fermions come from a quan- 
tum order in our vacuum. From the connection between 
string-net condensation, quantum order, and massless 
gauge/fermion excitations, it is very tempting to pro- 
pose the following answers to the fundamental questions 
about light and (nearly) massless fermions: 



What are light and fermions? 
Light is a fluctuation of condensed string-nets of arbi- 
trary sizes. Fermions are ends of open strings. 



Where light and (necirly) massless fermions 
come from? 

Light and the fermions come from the collective mo- 
tions of nets of string-like objects that fill our vacuum. 



Why light and (nearly) massless fermions exist? 
Light and the fermions exist because our vacuum 
chooses to have a string-net condensation. 



Had our vacuum chosen to have a "particle" conden- 
sation, thc;re would be only Nambu-Goldstone bosons 
at low energies. Such a universe would be very bor- 
ing. String-net condensation and the resulting light and 
(nearly) massless fermions provide a much more inter- 
esting universe, at least interesting enough to support 
intelligent life to study the origin of light and massless 
fermions. 

The string-net picture of fermions explains why there 
is always an even number of fermions in our universe. 
The string-net picture for gauge bosons and fermions also 
has an experimental prediction: all fermions must carry 
certain gauge charges. [14] At first sight, this prediction 
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appears to contradict with the known experimental fact 
that neutron carry no gauge charges. Thus one may think 
the string-net picture of gauge bosons and fermions has 
already been falsified by experiments. Here we would like 
to point out that the string-net picture of gauge bosons 
and fermions can still be correct if we assume the exis- 
tence of a new discrete gauge field, such as a gauge 
field, in our universe. In this case, neutrons and neutrinos 
carry a non-zero charge of the discrete gauge field. There- 
fore, the string- net picture of gauge bosons and fermions 
predict the existence of discrete gauge excitations (such 
as gauge flux lines) in our universe. 

We woiild like to remark that, despite the similarity, 
the above string-net picture of gauge bosons and fermions 
is different from the picture of standard superstring the- 
ory. In standard superstring theory, closed strings corre- 
spond to gravitons, and open string correspond to gauge 
bosons. All the elementary particles correspond to dif- 
ferent vibration modes of small strings in the superstring 
theory. Also, the fermions in the standard superstring 
theory come from the fermion fields on the world sheet. 
In our string-net picture, the vacuum is filled with large 
nets of strings. The massless gauge bosons correspond 
to the fiuctuations of large closed string-nets {ie nets of 
closed strings) and fermions correspond to the ends of 
open strings in string nets. Anti-commuting fields are 
not needed to produce (nearly) massless fermions. Mass- 
less fermions appear as low energy collective modes in a 
purely bosonic system. 

The string-net picture for gauge theories have a long 
history. The closed-string description of gauge fluctua- 
tions is intimately related to the Wilson loop in gauge 
theory. [17-19] The relation between dynamical gauge 
theory and a dynamical Wilson-loop theory was sug- 
gested in Ref. [20, 21]. Ref. [22] studied the Hamilto- 
nian of a non-local model - lattice gauge theory. It was 
found that the lattice gauge theory contains a string- 
net structure and the gauge charges can be viewed as 
ends of strings. In Ref. [23, 24] various duality relations 
between lattice gauge theories and theories of extended 
objects were reviewed. In particular, some statistical 
lattice gauge models were found to be dual to certain 
statistical membrane models. [25] This duality relation is 
directly connected to the relation between gauge theory 
and closed-string-net theory[13] in quantum models. 

Emerging fermions from local bosonic models also have 
a complicated history. The first examples of emerg- 
ing fermions/anyons were the fractional quantum Hall 
states, [4, 5] where fermionic/anyonic excitations were ob- 
tained theoretically from interacting bosons in magnetic 
field. [26] In 1987, fermion fields and gauge fields were 
introduced to express the spin- 1/2 Hamiltonian in the 
slave-boson approach. [27, 28] However, writing a bosonic 
Hamiltonian in terms of fermion fields does not imply 
the appearance of well defined fermionic quasiparticles. 
Emerging fermionic excitations can appear only in decon- 
fined phases of the gauge field. Ref. [29-32] constructed 
several deconfined phases where the fermion fields do de- 
scribe well defined quasiparticles. However, depending 
on the property of deconfined phases, those quasiparti- 
cles may carry fractional statistics (for the chiral spin 



states) [29, 30, 33] or Fermi statistics (for the Z2 decon- 
fined states). [31, 32] 

Also in 1987, in a study of resonating-valence-bond 
(RVB) states, emerging fermions (the spinous) were pro- 
posed in a nearest neighbor dimer model on square 
lattice. [34 36] But, according to the deconfinement pic- 
ture, the results in Ref. [34, 35] are valid only when the 
ground state of the dimer model is in the Z2 deconfined 
phase. It appears that the dimer liquid on square lattice 
with only nearest neighbor dimers is not a deconfined 
state, [35, 36] and thus it is not clear if the nearest neigh- 
bor dimer model on square lattice[35] has the fermionic 
quasiparticles or not. [36] However, on triangular lattice, 
the dimer liquid is indeed a Z2 deconfined state. [37] 
Therefore, the results in Ref. [34, 35] are valid for the 
triangular-lattice dimer model and fermionic quasiparti- 
cles do emerge in a dimer liquid on triangular lattice. 

All the above models with emerging fermions are 2-l-lD 
models, where the emerging fermions can be understood 
from binding flux to a charged particle. [26] Recently, it 
was pointed out in Ref. [14] that the key to emerging 
fermions is a string structure. Fermions can generally 
appear as ends of open strings. The string picture allows 
a construction of a 3-l-lD local bosonic model that has 
emerging fermions. 

Comparing with those previous results, the new fea- 
tures discussed in this paper are: (A) Massless gauge 
bosons and fermions can emerge from local bosonic mod- 
els as a result of string-net condensation. (B) Mass- 
less fermions are protected by the string-net condensa- 
tion (and the associated PSG). (C) String-net condensed 
states represent a new kind of phases which cannot be 
described Landau's symmetry breaking theory. Diff'erent 
string-net condensed states arc characterized by differ- 
ent PSG's. (D) QED and QCD can emerge from a local 
bosonic model on cubic lattice. The effective QED and 
QCD has AN families of leptons and quarks. Each family 
has one lepton and two flavors of quarks. 

The bottom line is that, within local bosonic mod- 
els, massless fermions do not just emerge by themselves. 
Emerging massless fermions, emerging massless gauge 
bosons, string-net condensations, and PSG are intimately 
related. They are just different sides of same coin - quan- 
tum order. 

According to the picture of quantum order, elemen- 
tary particles (such as photon and electron) may not be 
elementary after all. They may be collective excitations 
of a local bosonic system below Planck scale. Since we 
cannot do experiments close to Planck scale, it is hard 
to determine if photon and electron are elementary par- 
ticles or not. In this paper, we would like to show that 
the string-net picture of light and fermions is at least self 
consistent by studying some concrete local boson mod- 
els which produce massless gauge bosons and massless 
fermions through string-net condensations. The local bo- 
son models studied here are just a few examples among a 
long list of local boson modeis[8, 28, 29, 31-33, 35, 37-47] 
that contain emerging fermions and gauge fields. 

Here we would like to stress that the string-net pic- 
ture for the actual gauge bosons and fermions in our 
universe is only a proposal at moment. Although string- 
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net condensation can produce and protect massless pho- 
tons, gluons, quarks, and other charged leptons, we do 
not know at moment if string-net condensations can pro- 
duce neutrinos which are chiral fermions, and the weak- 
interaction SU{2) gauge field which couples chirally to 
the quarks and the leptons. Also, we do not know if 
string-net condensation can produce an odd number of 
families of quarks and leptons. The QED and QCD pro- 
duced by the known string-net condensations all contain 
an even number of families so far. The correctness of 
string-net condensation in our vacuum depend on resolv- 
ing the above problems. Nature has four fascinating and 
somewhat strange properties: gauge bosons, Fermi statis- 
tics, chiral fermions, and gravity. The string-net conden- 
sation picture provides a natural explanation for the first 
two properties. Two more to go. 

On the other hand, if we concern about a condensed 
matter problem: How to use bosons to make artificial 
light and artificial fermions, then the string- net picture 
and quantum order do provide an answer. To make ar- 
tificial light and artificial fermions, we simply let certain 
string-nets to condense. 

In some recent work, quantum orders and their connec- 
tion to emerging gauge bosons and fermions were studied 
using PSG's, without realizing their connection to string- 
net condensation. [9, 15, 46] In this paper, we will show 
that the quantum ordered states described by PSG's are 
actually string-net condensed states. The gauge bosons 
and fermions produced and protected by the PSG's have 
a very natural string-net interpretation. [13, 14] Quan- 
tum order, PSG, and string-net condensation are differ- 
ent parts of the same story. Here we will summarize and 
expand those previous work and try to present a coherent 
picture for quantum order, PSG, and string-net conden- 
sation, as well as the associated emerging gauge bosons 
and fermions. 



F. Organization 

Section III reviews the work in Ref. [14]. We will study 
an exactly soluble spin-1/2 model on square lattice[8, 46]. 
The model was solved using slave-boson approach. [46] 
This allowed us to identify the PSG that characterizes 
the non-trivial quantum order in the groimd state. [46] 
Here, following Ref. [14], we will solve the model from 
string-net condensation point of view. Since the ground 
state of the model can be described by both string-net 
condensation and PSG, this allows us to demonstrate 
the direct connection between string-net condensation 
and PSG in section IV. The model is also one of the 
simplest models that demonstrates the connection be- 
tween string-net condensation and emerging gauge field 
and fermions. [8, 14] 

However, the above exact soluble model does not con- 
tain gapless gauge boson and gapless fermions. If we 
regard the lattice scale as the Planck scale, then gauge 
bosons and fermions do not "exist" in our model in the 
sense discussed in section I A. In section V, we will dis- 
cuss an exact soluble local bosonic model that contain 
massless Dirac fermions. In sections VII and VIII, we 



will discuss local bosonic models that give rise to massless 
electrons, quarks, gluons, and photons. Gauge bosons 
and fermions "exist" in those latter models. 



II. LOCAL BOSONIC MODELS 

In this paper, we will only consider local bosonic mod- 
els. Local bosonic models are important since they are 
really local. Wc note that a fermionic model arc in gen- 
eral non local since the fermion operators at different sites 
do not commute, even when the sites are well separated. 
Due to their intrinsic locality, local bosonic models are 
natural candidates for the fundamental theory of nature. 
In the following we will give a detailed definition of local 
bosonic models. 

To define a physical system, wc need to specify (A) 
a total Hilbert space, (B) a dcfiuition of a set of local 
physical operators, and (C) a Hamiltonian. With this 
understanding, a local bosonic model is defined to be a 
model that satisfies: (A) The total Hilbert space is a di- 
rect product of local Hilbert spaces of finite dimensions. 
(B) Local physical operators are local bosonic operators. 
By definition, local bosonic operators arc operators acting 
within a local Hilbert space or finite products of those op- 
erators for nearby local Hilbert spaces. Those operators 
are called local bosonic operators since they all commute 
with each other when far apart. (C) The Hamiltonian is 
a sum of local physical operators. 

A spin-1/2 system on a lattice is an example of local 
bosonic models. The local Hilbert space is two dimen- 
sional which contains | f) and | J,) states. Local physical 
operators are cr?, (yi'y'l_^_^, etc , where a"" , a = x,y,z are 
the Pauli matrices. 

A free spinless fermion system (in 2 or higher dimen- 
sions) is not a local bosonic model despite it has the 
same total Hilbert space as the spin-1/2 system. This 
is because the fermion operators Cj on different sites do 
not commute and are not local bosonic operators. More 
importantly, the fermion hoping Hamiltonian in 2 and 
higher dimensions cannot be written as a sum of local 
bosonic operators. (Note in higher dimensions, we can- 
not write all the hoping terms ctcj as product of local 
bosonic operators. However, due to the Jordan- Wigner 
transformation, a ID fermion hoping c|_|_]^Ci can be writ- 
ten as a local bosonic operators. Hence, a ID fermion 
system can be a local bosonic model if we exclude Ci 
from our definition of local physical operators.) 

The bosonic field theory without cut-off is not a local 
bosonic model. This is because the local Hilbert space 
does not have a finite dimension. A lattice gauge theory 
is not a local bosonic model. This is because its total 
Hilbert space cannot be a direct product of local Hilbert 
spaces. 

Another counter example of local bosonic model is a 
quantum closed-string-net model. A quantum closed- 
string-net model on lattice can be defined in the follow- 
ing way. Let us consider only strings that cover nearest 
neighbor links. A closed-string configuration may have 
many closed strings with or without overlap. We will 
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even plaquette 



FIG. 2: An open-string excitation on top of the ground state 
of Hj. 

caU a closed-string configuration a closed string- net. For 
every closed string-net , we assign a quantum state. All 
those quantum states form a basis of the total Hilbert 
space of the closed-string-net model. Just like lattice 
gauge theory, the closed-string-net model is not a local 
bosonic model since the total Hilbert space cannot be 
a direct product of local Hilbert spaces. It turns out 
that closed-string-net models and lattice gauge model 
are closely related. In fact some closed-string-net models 
(or statistical membrane models) are equivalent to lattice 
gauge models. [22-25] 

III. Z2 SPIN LIQUID AND STRING-NET 
CONDENSATION ON SQUARE LATTICE 

A. Hamiltonians with closed-string-net 
condensation 

Let us first consider an arbitrary spin-1/2 model on a 
square lattice. The first question that we want to ask 
is what kind of spin interaction can given rise to a low 
energy gauge theory. If we believe the connection be- 
tween gauge theory and closed-string-net theory, [13, 22- 
25] then one way to obtain a low energy gauge theory 
is to design a spin interaction that allow strong fluctua- 
tions of large closed string-nets, but forbid other types of 
fluctuations (such as local spin flip, open string-net fluc- 
tuations, etc ). (Note that closed string-nets are nets of 
strings formed by intersecting/overlapping closed strings, 
while open string-nets are nets of strings containing at 
least an open string.) We hope the presence of strong 
fluctuations of large closed-strings will lead to conden- 
sation of closed strings of arbitrary sizes, which in turn 
gives rise to a low energy gauge theory. 

Let us start with 

Hj = -jY.crt-jY,ay (1) 

even odd 

where i = {ix,iy) labels the lattice sites, cr^'^'^ are the 
Pauli matrices, and J2even Sorfd) ^ '-'^^^ even 
sites with (— )* = (— 1)*^+'« — 1 (or over odd sites with 
{-Y EE {-iy-+iy = -1). The ground state of Hj, |0), 
has spins pointing to x-direction on even sites and to y- 
direction on odd sites (see Fig. 2). Such a state will be 
defined as a state with no string. 



To create a string excitation, we first draw a string 
that connect nearest neighbor even plaquettes (see Fig. 
2). We then flip the spins in the string. Such a string 
state is created by the following string creation operator 
(or simply, string operator): 

WiC) = l[cTr (2) 
c 

where the product Yic is over all the sites on the string, 
tti = y if i is even and ai — x if i is odd. A generic string 
state has a form 

\CiC2...) = W{Ci)W{C2)...\Q) (3) 

where Ci, C'2, ■■■ are strings with no overlapping ends. 
Such a state will be called a string-net state and 

W{Cnet) = W{Ci)W{C2)... 

will be called a string-net operator. The state IC1C2...) 
is an open-string-net state if at least one of Ci is an 
open string. The corresponding operator W{Cnet) will 
be called an open-string-net operator. If all Ci are closed 
loops, then IC1C2...) is an closed-string-net state and 
W{Cnet) an closed-string-net operator. The Hamiltonian 
has no string-net condensation since its ground state |0) 
contains no string-nets. To obtain a Hamiltonian with 
closed-string-net condensation, we need to first find a 
Hamiltonian whose ground state contains a lot of closed 
string-nets of arbitrary sizes and do not contain open 
string-nets. 

Let us first write down a Hamiltonian such that closed 
strings cost no energy and any open strings cost a large 
energy. One such Hamiltonian has a form 

even 

Fi =<<yl+^<+.+y<yt+y (4) 

We find the no-string state |0) is one of the ground state 
of Hu (assuming [/ > 0) with energy —UNsHe- AH the 
closed-string-net states, such as W{Cciose)\Q), are also 
ground state of Hu since [Hu ,W{Cciose)\ = 0. An open- 
string state W(Copen )|0) is also an eigenstate of Hu but 
with energy —UNsHe + 2C/. We see that each end of open 
string cost an energy U. We also note that the energy 
of closed strings does not depend on the length of closed 
strings. Thus the closed strings in Hu have no tension. 
We can introduce a string tension by adding the Hj to 
our Hamiltonian. The string tension will be 2J per site 
(or per segment). We note that, any string-net state 
IC1C2...) is an eigenstate of Hu + Hj. Thus, string-nets 
in the model described by Hu + Hj do not fluctuate and 
hence cannot condense. To make string-nets to fluctuate, 
we need a g-term 

Hg^gJ2u{Cp) (5) 
p 

where p labels the odd plaquettes and Cp is the closed 
string around the plaquette p. In fact 

Hg = -gY,Fi (6) 

odd 
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FIG. 3: The proposed phase diagram for the H = Hjj + 
Hg + Hj model. J is assumed to be positive. The four 
string-net condensed phases are characterized by a pair of 
PSG's {PSGcharge,PSGvortex)- MO marks au magnetic or- 
dered state. 

This way, we obtain the Hamiltonian of our spin-1/2 
model 

H^Hu + Hj + Hg (7) 

B. String condensation and low energy effective 
theory 

When J = in Eq. (7), the model is exactly soluble 
since [Fj, Fj] 0.[8, 46] All the eigenstates of Hu + Hg 
can be obtained from the common eigenstates of Fi. 
Since F? = 1, the eigenvalues of Fi are simply ±1. Thu 
s all the eigenstates of Hu + Hg are labeled by ±1 on 
each plaquette. (Note, this is not true for finite sys- 
tems where the boundary condition introduce additional 
complications. [46]) The energies of those eigenstates are 
sum of eigenvalues of Fi weighted by U and g. 

From the results of exact soluble model, we suggest a 
phase diagram of our model as sketched in Fig. 3. We 
will show that the phase diagram contains four different 
string-net condensed phases and one phase with no string 
condensation. All the phases have the same symmetry 
and are distinguished only by their different quantum 
orders. 

Let us first discuss the phase with U,g > 0. We will 
assume J = and U ^ g. In this limit, all states con- 
taining open strings will have an energy of order U . The 
low energy states contain only closed strings (or more 
generally closed string-nets) and satisfy 

Fi\i—even — 1 (8) 

For infinite systems, the different low energy states are 
labeled by the eigenvalues of Fi on odd plaquettes: 

Fi\i^odd^±l (9) 

In particular, the ground state is given by 

FiU^odd = 1. (10) 
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FIG. 4: A hopping of the Z2 charge around four nearest neigh- 
bor even plaquettes. 

All the closed-string-net operators W{Cnet) commute 
with Hjj + Hg. Hence the ground state \'^o) of Hu + Hg 
satisfies 

(*0|W^(C„et)|*0) = 1. (11) 

Thus the U,g > ground state has a closed-stringi- 
net condensation. The low energy excitations above the 
ground state can be obtained by flipping Fi from 1 to —1 
on some odd plaquettes. 

If we view Fi on odd plaquettes as the flux in Z2 gauge 
theory, we find that the low energy sector of model is 
identical to a Z2 lattice gauge theory, at least for infi- 
nite systems. This suggests that the low energy effective 
theory of our model is a Z2 lattice gauge theory. 

However, one may object this result by pointing out 
that the low energy sector of our model is also identical to 
an Ising model with one spin on each the odd plaquette. 
Thus the the low energy effective theory should be the 
Ising model. We would like to point out that although 
the low energy sector of our model is identical to an Ising 
model for infinite systems, the low energy sector of our 
model is different from an Ising model for finite systems. 
For example, on a finite even by even lattice with periodic 
boundary condition, the ground state of our model has 
a four- fold degeneracy. [8, 46] The Ising model does not 
have such a degeneracy. Also, our model contains an 
excitation that can be identified as Z2 charge (see below) . 
Therefore, the low energy effective theory of our model is 
a Z2 lattice gauge theory instead of an Ising model. The 
Fi — —1 excitations on odd plaquettes can be viewed as 
the Z2 vortex excitations in the Z2 lattice gauge theory. 

C. Three types of strings and emerging fermions 

What is the Z2 charge excitations? We note that, in 
the closed-string-net condensed state, the action of the 
closed-string operator Eq. (2) on the ground state is triv- 
ial. This suggests that the action of the open-string op- 
erators on the ground state only depend on the ends of 
strings, since two open strings with the same ends only 
differ by a closed string. Therefore, an open-string op- 
erator create two particles at its ends when acting on 
the string condensed state. Since the strings in Eq. (2) 
only connect even plaquettes, the particle corresponding 
to the ends of the open strings always live on the even 
plaquettes. We will call such a string Tl string. Form 
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the commutation relation between Fi and the open-string 
operators, we find that the open-string operators flip the 
sign of Fi at its ends. Thus each particle created by 
the open-string operators has an energy 2U. Now, let 
us consider the hopping of one such particle around four 
nearest neighbor even plaquettes (see Fig. 4). We see 
that the product of the the four hopping amplitudes is 
given by the eigenvalue of Fi on the odd plaquette in the 
middle of the four even plaquettes. [8, 14] This is exactly 
the relation between charge and flux. Thus if we identify 
Fi on odd plaquettes as Z2 flux, then the ends of strings 
on even plaquettes will correspond to the Z2 charges. 
We note that, due to the closed-string condensation, the 
ends of open strings are not confined and have only short 
ranged interactions between them. Thus the Z2 charges 
behave like quasiparticles with no string attached. 

Just like the Z2 charges, a pair Z2 vortices is also cre- 
ated by an open string operator. Since the Z2 vortices 
correspond to flipped Fi on odd plaquettes, the open- 
string operator that create Z2 vortices is also given by 
Eq. (2), except now the product is over a string that con- 
nect odd plaquettes. We will call such a string T2 string. 
(The strings connecting even plaquettes were called Tl 
strings.) 

We would like to point out that the reference state 
{ie the no string state) for the T2 string is different from 
that of the Tl string. The no-T2-string state is given 
by |0) with spin pointing in y-direction on even sites and 
x-direction on odd sites. Since the Tl and T2 strings 
have different reference state, we cannot have a dilute gas 
of the Tl strings and the T2 strings at the same time. 
One can easily check that the T2 string operators also 
commute with Hu + Hg. Therefore, the ground state 
\^o)i in addition to the Tl closed-string condensation, 
also has a T2 closed-string condensation. 

The hopping of a Z2 vortex is induced by a short T2 
open-string. Since the T2 open-strings operators all com- 
mute with each other, the Z2 vortex behave like bosons. 
Similarly, the Z2 charges also behave like bosons. How- 
ever, Tl open-string operators and T2 open-string oper- 
ators do not commute. As a result, the ends of Tl string 
and the ends of T2 string have nontrivial mutual statis- 
tics. As we have already shown that moving a Z2 charge 
around a Z2 vortex generate a phase tt, the Z2 charges 
and the Z2 vortices have a semionic mutual statistics. 

The T3 strings are defined as bound states of Tl and 
T2 strings. The T3 string operator has a form W(C) = 
Yin '^iZ ■ where C is a string connecting the mid-points of 
the neighboring links (see Fig. 6), and i„ are sites on the 
string. Ijn = z if the string does not turn at site im (see 
Fig. 6). /m = X or y if the string makes a turn at site 
im- Im = X ii the turn forms a upper-right or lower- left 
corner. = y if the turn forms a lower-right or upper- 
left corner. (See Fig. 6.) The ground state also has a 
condensation of T3 closed-strings. The ends of T3 string, 
as bound states of the Z2 charges and the Z2 vortices, 
are fermions. The bound state is formed by a Z2 charge 
and a Z2 vortex on the two plaquettes on the two side 
of a link (ie Fi = —1 on the two sides of the hnk). Thus 
the fermions live on the links. It is interesting to see that 



string-net condensation in our model directly leads to Z2 
gauge structure and three new type of quasiparticles: Z2 
charge, Z2 vortex, and fermions. Fermions, as ends of 
open T3 strings, emerge from our purely bosonic model. 

Since ends of Tl string are Z2 charges, the Tl string 
can be viewed as strings of Z2 "electric" flux. Similarly, 
the T2 string can be viewed as strings of Z2 "magnetic" 
flux. 



IV. CLASSIFICATION OF DIFFERENT 
STRING CONDENSATIONS BY PSG 

A. Four classes of string-net condensations 

As we have seen in last section that when U > 0, g > 0, 
and J = 0, the ground state of our model is given by 

Pi\i=even — Ij Pi\i=odd — !• (l^) 

We will call such a phase Z2 phase to stress the low en- 
ergy Z2 gauge structure. In the Z2 phase, the Tl string 
operator Wx{C\) and the T2 string operator W2{C2) have 
the following expectation values 

(W^i(Ci)) = 1, (1^2(^2)) = 1 (13) 

When U > 0, g < 0, and J = 0, the ground state is 
given by 

Fi\i—even — 1, Fi\i—Q(i(i — 1. (f^) 

We see that there is 7r-flux through each odd plaquette. 
We will call such a phase Z2-flux phase. The Tl string 
operator and the T2 string operator have the following 
expectation values 

(T4^i(Ci)) = (-)^°-, {W2{C2)) = 1 (15) 

where Nodd is the number of odd-plaquettes enclosed by 
the Tl string Ci. 

When U < 0, 5 > 0, and J = 0, the ground state is 

Pi\i=even — 1; Pi\i=odd — !• (i-^) 

The ground state has a Z2 charge on each even plaquette. 
We will call such a phase Z2-charge phase. The Tl string 
operator and the T2 string operator have the following 
expectation values 

(VFi(Ci)) = 1, {W2{C2)) = (-)^— (17) 

where Neven is the number of even-plaquettes enclosed by 
the T2 string C2. Note that the Z2-flux phase and the 
Z2-charge phase, different only by a lattice translation, 
are essentially the same phase. 

When U < 0, g < 0, and J = 0, the ground state 
becomes 

Pi\i=even — 1; Pi\i=odd — !• (1^) 

There is a Z2 charge on each even plaquette; and 7r-flux 
through each odd plaquette. We will call such a phase 
Z2-flux-charge phase. The Tl string operator and the T2 
string operator have the following expectation values 

(Wi(Ci)) = (-)^"<^^ (W2(C2)) = (-)^— (19) 
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B. PSG and ends of condensed strings 

Prom the different {Wi (Ci )) and ( VF2 (C2 )) , we see that 
the above four phases have different string-net condensa- 
tions. However, they ah have the same symmetry. This 
raises an issue. Without symmetry breaking, how do we 
know the above four phases are really different phases? 
How do we know that it is impossible to change one 
string-net condensed state to another without a phase 
transition? 

In the following, we will show that the different string- 
net condensations can be described by different PSG's 
(just like different symmetry breaking orders can be de- 
scribed by different symmetry groups of ground states.) 
In Ref. [9; 10], different quantmn orders were introduced 
via their different PSG's. The connection between string- 
net condensation and PSG allows us to connect string- 
net condensation to the quantum order introduced in 
Ref. [9, 10]. In particular, the PSG's are shown to be 
a universal property of a quantum phase, which can be 
changed only by phase transitions. Thus the different 
PSG's for the different string-net condensed states indi- 
cate that those different string-net condensed states be- 
long to different quantum phases. 

When closed-string-nets condense, the ends of open 
strings behave like independent particles. Let us consider 
two particles states \P1P2) described by the two ends of a 
Tl string. Note that the ends of the Tl strings, and hence 
the Z2 charges, only live on the even plaquettes. Here Pi 
and P2 label the even plaquettes. For our model Hjj+Hg, 
\P1P2) is an energy eigenstate and the Z2 charges do not 
hop. Here we would like to add a term 



(20) 



to the Hamiltonian. The i-term t^^(<j^ + af) makes 
the Z2 charges to hop among the even plaquettes with a 
hopping amplitude of order t. The dynamics of the two 
Z2 charges is described by the following effective Hamil- 
tonian in the two-particle Hilbert space 



H = H{p,) + Hip2) 



(21) 



where H{pi) describes the hopping of the first particle Pi 
and -ff (P2) describes the hopping of the second particle 
P2- Now we can define the PSG in a string-net condensed 
state. The PSG is nothing but the symmetry group of 
the hopping Hamiltonian H{p). 

Due to the translation symmetry of the underlying 
model Hu + Hg + Ht, we may naively expect the hopping 
Hamiltonian of the Z2 charge H{p) also have a transla- 
tion symmetry 

H{P) =TlH{p)T,y, T,y\p) = \p + x + y) 



H{p) =TlyH{p)T, 



xyl 

Txy\p) = \p + x-y) (22) 



The above implies PSG = translation symmetry group. 

It turns out that Eq. (22) is too strong. The underlying 
spin model can have translation symmetry even when 
H{p) does not satisfy Eq. (22). However, the possible 



symmetry groups of H{p) (the PSG's) are strongly con- 
strained by the translation symmetry of the underlying 
spin model. In the follow, we will explain why the PSG 
can be different from the symmetry group of the physi- 
cal spin model, and what conditions that the PSG must 
satisfy in order to be consistent with the translation sym- 
metry of the spin model. 

We note that a string always has two ends. Thus a 
physical state always has an even number of Z2 charges. 
The actions of translation on a two-particle state are 
given by 



Ti^^ \Pi , P2) =e^«''(Pi'P^) \p,+x + y,P2 + x + y) 

y) (23) 



Ti^^\Pi,P2) =e^==--(f-P=^)|pi +x-y,p2 + x 



The phases e^^^f^^i'^^) and e^='s(Pi.P2) come from the am- 
biguity of the location of the string that connect p^ and 

P2. ie the phases can be different if the string connecting 

(2) (2) 

the two Z2 charges has different locations. Txy and T^y' 
satisfy the algebra of translations 



7^(2)7.(2) 

■^xy ^xy 



. rp{2)rp(2) 

-^xy -^xy 



(24) 



Txy"^ and T^^' are direct products of translation opera- 
tors on the single-particle states. Thus, in some sense, 
the single-particle translations are square roots of two- 
particle translations. 

The most general form of single-particle translations 
is given by T^yGxy and T^yGxy, where the actions of 
operators T^y^xy and Gxy,xy are defined as 

T,y\p) =\p + x + y) 
Txy\p) =\p + x-y) 
Gxy\p) =e^^-(P)|p) 



.(2) 



Gxy\p} 



p) 



(25) 



In order for the direct product ri^-* = T^yGxy ® T^yG, 



(2) 



xy 

TxyGxy®TxyGxy to rcproducc the translation 
algebra Eq. (24), we only require T^yGxy and T^yGxy to 



find T^y 



satisfy 



or 



^xyG xy^xyG xy — ^xy^xy^xy^xy 



'^xyG xyTxyG xy T-j-yG xyTxyG xy 



(26) 



(27) 



The operators T^yGxy and T^yGx-g generate a group. 
Such a group is the PSG introduced in Ref. [9]. The 
two different algebra Eq. (26) and Eq. (27) generate two 
different PSG's, both are consistent with the translation 
group acting on the two-particle states. We will call the 
PSG generated by Eq. (26) Z2A PSG and the PSG gen- 
erated by Eq. (27) Z2B PSG. 

Let us give a more general definition of PSG. A PSG is 
a group. It is a extension of symmetry group (SG), ie a 
PSG contain a normal subgroup (called invariant gauge 
group or IGG) such that 



PSG/IGG = SG 



(28) 
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For our case, the SG is the translation group SG = 
{l,T^l\Tif ,...}. For every element in SG, a^^) e SG, 
there are one or several elements in PSG, a G PSG, such 
that a® a = a^'^\ The IGG in our PSG is formed by 
the transformations Go on the singe-particle states that 
satisfy Go (S) Go = 1- We find that IGG is generated by 

Go\p) = -\p> (29) 

Go, TxyGxy and T^yGxy generate the Z2A and Z2B 
PSG's. 

Now we see that the underlying translation symmetry 
does not require the single-particle hopping Hamiltonian 
H{p) to have a translation symmetry. It only require 
h{p) to be invariant under the Z2A PSG or the Z2B 
PSG. When H{p) is invariant under the Z2A PSG, the 
hopping Hamiltonian has the usual translation symme- 
try. When H(j)) is invariant under the Z2B PSG, the 
hopping Hamiltonian has a magnetic translation symme- 
try describing a hopping in a magnetic field with tt-Aux 
through each odd plaquette. 

C. PSG's classify different string-net condensations 

After understand the possible PSG's for the hopping 
Hamiltonian of the ends of strings, now we are ready to 

calculate the actual PSG's. Let us consider two ground 
states of our model Hu + Hg+Ht. One has -Fi |i=odd = 1 
(for g > 0) and the other has Fi\i^odd = — 1 (for g < 0). 
Both ground states have the same translation symmetry 
in x+y and x—y directions. However, the corresponding 
single-particle hopping Hamiltonian H{p) has different 
symmetries. For the Fi\i=odd = 1 state, there is no flux 
through odd plaquettes and H{p) has the usual transla- 
tion symmetry. It is invariant under the Z2A PSG. While 
for the Fi\i=odd = ^1 state, there is 7r-flux through odd 
plaquettes and H{p) has a magnetic translation symme- 
try. Its PSG is the Z2B PSG. Thus the Fi\i=odd = 1 state 
and the Fi\i=odd = ^1 state have different orders despite 
they have the same symmetry. The different quantum 
orders in the two states can be characterized by their 
different PSG's. 

The above discussion also apply to the Z2 vortex and 
T2 strings. Thus the quantum orders in our model are 
described by a pair of PSG's (PSGcharge, PSGyoriex), 
one for the Z2 charge and one for the Z2 vortex. The 

PSG pairs {PSGchargei P SGyortex 

) allows us to distin- 
guish four different string-net condensed states of model 
H = Hu + Hg + Ht.{See Fig. 5.) 

Now let us assume U = g in our model: 

Hu + Hg + Ht = Ht-VY,Fi (30) 

i 

The new physical spin model has a larger translation 
symmetry generated by Ai = a; and Ai = y (see Fig. 
5). Due to the enlarged symmetry group, the quantum 
orders in the new system should be characterized by a 
new PSG. In the following, we will calculate the new 
PSG. 
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FIG. 5: The proposed phase diagram for the H = Hu + Hg + 
Ht model, t = t' is assumed to be positive. The four string- 
net condensed phases arc characterized by a pair of PSG's 
{PSGcharge,PSGvortex)- FM marlcs a ferromagnetic phase. 

The single-particle states are given by \p) . When p is 

even, |p) corresponds to a Z2 charge and when p is odd, 
\p) corresponds to a Z2 vortex. We see that a transla- 
tion by X (or y) will change a Z2 charge to a Z2 vortex 
or a Z2 vortex to a Z2 charge. Therefore the effective 
single-particle hopping Hamiltonian H{p) only contain 
hops between even plaquettes or odd plaquettes. The 
single-particle Hamiltonian H{p) is invariant under the 
following two transformations Go and Gq: 

Go|p) = -|p>, G'^\p) = (-r\p> (31) 

We note that Gq®Gq = G'q®G'q = I. Therefore both Go 
and Gq correspond to the identity element of the sym- 
metry group of two-particle states. (Go, Gq) generate the 
IGG of the new PSG. The new IGG is Z2 x Z2. 

The translations of single-particle states by x and by 
y are generated by T^Gx and TyGy. The translation by 
x + y and hy x — y are given by 

PxyGxy — PyGyTxGx 

FxyGxy ~{PyGy) TxGx (^^2) 

Since TxyGxy and TxyGxy are the translations of the Z2 
charge and the Z2 vortex discussed above, we find 

(F'xyGxy^ {'Fxy^xy^ 'FxyG xyTxyG xy — TJ (^3) 

where t? = 1 for the {Z2A, Z2A) state with Fi = 1 and 
T/ = -1 for the {Z2B,Z2B) state with Fi = -1. Also 

TxGx and TyGy must satisfy 

{TyGy)-\TxGx)-%GyTxGx e IGG (34) 
since on the two-particle states 

(Tf)-i(T(2))-^Tfr(2) = l (35) 

Therefore, {TyGy)-HTxGx)^^TyGyTxGx may take the 
following possible values 1, —1, (— )^, and — (— )p. Only 
choices rj^ and —rj^ are consistent with Eq. (33) and we 
have 

{TyGy)-\TxGx)-^TyGyTxGx = v' (36) 
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FIG. 6: Fermion hopping around a plaquette, around a 
square, and around a site. 



We like to point out that the different choices of 77' = 
±1 do not lead to different PSG's. This is because if 
TxGx is a symmetry of the H{p), then TxGx{—)^ is also 
a symmetry of the H{p). However, the change Gx — > 
Gx{—)^ will change the sign of 77'. Thus i]' — 1 and 
77' = —1 will lead to the same PSG. But the different 
signs of 77 will lead to different PSG's. 

iGo,G'o) and {TxGx,TyGy) generate the new PSG. The 
single-particle Hamiltonian H{p) is invariant under such 
a PSG. 77 = 1 and 77 = —1 correspond to two different 
PSG's that characterize two different quantum orders. 
The ground state for V >Q and \V\ > t (see Eq. (30)) is 
described by the 77 = 1 PSG. The ground state for V <0 
and \V\ > t is described by the 7; = -1 PSG. The two 
ground states have different quantum orders and different 
string-net condensations. 



D. Different PSG's from the ends of different 
condensed strings 

In this section we still assume U — g and consider only 
the translation invariant model Eq. (30). In the above we 
discussed the PSG for the ends of one type of condensed 
strings in different states. In this section, we will concen- 
trate on only one ground state. We know that the ground 
state of our spin- 1/2 model contain condensations of sev- 
eral type of strings. We like to calculate the the different 
PSG's for the different condensed strings. 

The PSG's for the condensed Tl and T2 strings were 
obtained above. Here we will discuss the PSG for the 
T3 string. Since the ends of the T3 strings live on the 
links, the corresponding single-particle hopping Hamil- 
tonian Hf{l) describes fermion hopping between links. 
Clearly, the symmetry group (the PSG) of Hf{l) can be 
different from that of H{p). 

Let us consider fermion hopping around some small 
loops. The four hops of a fermion around a site i (see Fig. 
6) are generated by erf, af, a^, and af. The total ampli- 
tude of a fermion hopping around a site is <^\<^i'^^'^i = 
— 1. The fermion hopping around a site always sees tt- 
fiux. The four hops of a fermion around a plaquette p (see 
Fig. 6) are generated by cr? , crf^+^, and al^^y, 

where io is the lower left corner of the plaquette p. The 
total amplitude of a fermion hopping around a plaquette 
is given by = i^io- When y > 0, 

the ground state has Fi = \. However, since site io is 
next to the end of T3 string, we have Fi^ — —Fi = —1. 



In this case, the fermion hopping around a plaquette sees 
TT-flux. For V < ground state, we find that fermion 
hopping around a plaquette sees no flux. 

Let us define the fermion hopping I I + x as the 
combination of two hops I ^ I + ^ — ^ I + x and 
the fermion hopping I I + y as the combination of 
1^1 + ^ + ^ ^l + y (see Fig. 6). Under such a 
definition, a fermion hopping around a square I I + 
x~^l + x + y^l + y^l correspond to a fermion 
hopping around a site and a fermion hopping around a 
plaquette discussed above (see Fig. 6). Therefore, the 
total amplitude for a fermion hopping around a square is 
given by the sign of V: sgn{V). We find the translation 
symmetries {TxGx,TyGy) of the fermion hopping Hf(l) 
satisfies 

iTyGy)-\TxGx)-%GyTxGx = sgn(y) (37) 

which is different from the translation algebra for H{p) 
Eq. (36). Hf{l) is also invariant under Gq: 
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\l) 



(38) 



{GQ,TxGx,TyGy) generate the symmetry group - the 
fermion PSG - of Hf{l). We will call the fermion PSG 
Eq. (37) for sgn(y) = 1 the Z2A PSG and the fermion 
PSG for sgn(y) = -1 the Z2B PSG. We see that the 
quantum orders in the ground state can also be charac- 
terized using the fermion PSG. The quantum order in 
the V > Q ground state is characterized by the Z2A PSG 
and the quantum order in the y < ground state is 
characterized by the Z2B PSG. 

In Ref. [46], the spin-1/2 model Eq. (30) (with t = 
t' = 0) was viewed as a hardcore boson model. The 
model was solved using slave-boson approach by split- 
ting a boson into two fermions. Then it was shown the 
fermion hopping Hamiltonian for V > Q and V <0 states 
have different symmetries, or invariant under different 
PSG's. According to the arguments in Ref. [9], the dif- 
ferent PSG's imply different quantum orders in the V > 
and V < Q ground state states. The PSG's obtained in 
Ref. [46] for the ^ > and < phases agrees exactly 
with the fermion PSG's that we obtained above. This ex- 
ample shows that the PSG's introduced in Ref. [10, 46] 
are the symmetry groups of the hopping Hamiltonian of 
the ends of condensed strings. The PSG description and 
the string-net-condensation description of quantum or- 
ders are intimately related. 

Here we would like to point out that the PSG's intro- 
duced in Ref. [9, 10] are all fermion PSG's. They are only 
one of many different kinds of PSG's that can be used to 
characterize quantum orders. In general, a quantum or- 
dered state may contain condensations of several types of 
strings. The ends of each type of condensed strings will 
have their own PSG. 



V. MASSLESS FERMION AND PSG IN 
STRING-NET CONDENSED STATE 

In Ref. [9, 16], it was pointed out that PSG can pro- 
tect masslessness of the emerging fermions, just like sym- 
metry can protect the masslessness of Nambu-Goldstone 
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bosons. In this section, wc arc going to study an ex- 
act soluble spin-ii model with string- net condensation 
and emerging massless fermions. Through this soluble 
model, we demonstrate how PSG that characterizes the 
string-net condensation can protect the masslessness of 
the fermions. The exact soluble model that we are going 
to study is motivated by Kitaev's exact soluble spin- 1/2 
model on honeycomb lattice. [48] 



A. Exact soluble spin-ii model 

The exact soluble model is a local bosonic model on 
square lattice. To construct the model, we start with four 
Majorana fermions A?, a = x,x,y,y and one complex 
fermion -0. A? satisiy 



We note that 



Ui,i+y — -iXiXl 



Ui 



(39) 



(40) 



form a commuting set of operators. Using such a com 

muting set of operators, we can c;onstruct the following 
exact soluble interacting fermion model 

H =g^Fi+t '^{iUi,i+x-4'i'4'i+x + iUi,i+yil>lil>i+y + h.c. 

i i 

where ii = i + x, 12 = i + x + y, = i + and t is 
real. We will call Fi a Z2 flux operator. To obtain the 
Hilbert space within which the Hamiltonian H acts, we 
group A'^'^'^'^ into two complex fermion operators 



2^1 



X^i+iXl 2i;2,i = Xy+iXj 



(42) 



on each site. The complex fermion operators ^^1,2 and tp 
generate an eight dimensional Hilbert space on each site. 

Since Uij commute with each other, we can find the 
common eigenstates of the Uij operators: \{sij},n), 
where Sij is the eigenvalue of Uij, and n labels different 
degenerate common eigenstates . Since {Uij)"^ = 1 and 



Uij = Uji , Sij satisfies Sij = ±1 and Si 



Within 



the subspace with a fixed set of Sij-. {\{sij},n}\n = 
1,2,...}, the Hamiltonian has a form 

H =g^ fi + t ^(zsi,i+3,Vi^i+x + isi,i+yipiipi+y + h.c.) 

i i 

which is a free fermion Hamiltonian. Thus we can find all 
the many-body eigenstates of ipi: |{sij},^'„) and their 
energies E{{sij}, n) in each subspace. This way we solve 
the interacting fermion model exactly. 

We note that the Hamiltonian H can only change the 
fermion number on each site by an even number. Thus 
the H acts within a subspace which has an even num- 
ber of fermions on each site. We will call the subspace 
physical Hilbert space. The physical Hilbert space has 



only four states per site. When defined on the physical 
space, H becomes a local bosonic system which actually 
describes a spin \ y-\ system (with no spin rotation sym- 
metry). We will call such a system spin-ii system. To 
obtain an expression of H within the physical Hilbert 
space, we introduce two Majorana fermions ?7i j and 772,1 
to represent ipi: 2t/;i = r]i^i + iri2,i- We note that X°"qi, 
a = x,x,y,y, act within the four dimensional physical 
Hilbert space on each site, and thus are 4 by 4 matrices. 
Also {—iX°'r]i, —iX^'rji} = 2Sab, thus the four 4 by 4 matri- 
ces A"ryi satisfy the algebra of Dirac matrices. Therefore 
we can express A"r/i in terms of Dirac matrices 7": 



7^ =0-^ O C7^, 



7« =cr^ cr-^, 
We can also define the 7^ 



7^ =0-^ (g) (7^ 

7^ =a^ (g) (T^ 



=A^A*A^A^ = irjim 



(44) 



(45) 



where we have used 1 — 2tjj^ip = —ir]\r]2 and 
{—iX^X^){~iXyX^){—irjirj2) = 1 for states with even 
numbers of fermions. With the above definition of 7"^ 
and 7^, we find that 



a 5 



and 



X'^m = 7^7 



A>=-(7« + 7V) = n-^ 



(46) 



7-'" =(7+'")t 



We also have 



A" A'' = 7^7'' = i 



ah 



(47) 



(48) 



The above relations allows us to write H in terms of 4 
by 4 Dirac matrices. For example 



.^J/a: ay vx xy 
H H+x H+x+y H+y 



and 



*7i li+x 



Ui,i+xtlJitlJi+x ■■ 
Ui,i+yll^ii^i+y = - ilt'^^i+y 



(49) 



(50) 



The physical states in the physical Hilbert space are 
invariant under local Z2 gauge transformations generated 
by 



i 

rii =V'J_iVi,i -I- 4>\^i'4>2,i + V'lV'i 



(51) 



where Gi is an arbitrary function with only two values 
±1 and Hi the number of fermions on site i. We note 
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that the Zi gauge transformations change i/'/i ~^ Giil}ii. 
The projection into the physical Hilbcrt space with even 
fermion per site makes our theory a Z2 gauge theory. 

Since the Hamiltonian H in Eq. (41) is Z2 gauge in- 
variant: [G, i?] — 0, cigcnstate of H within the phys- 
ical Hilbert space can be obtained from |{sij},^'n) by 
projecting into the physical Hilbert space: 7'|{sij}, 
The projected state V\{sij], ^'„) (or the physical state), 
if non-zero, is an eigenstate of the spin-i^ model with en- 
ergy i?({sjj}, n). The Z2 gauge invariance implies that 



if Sij and s^j are Z2 gauge equivalent 



(52) 



(53) 



Let us count the states to show that the projected 
states 'P|{sij},^n) generate all states in the physical 
Hilbert space. Let us consider a periodic lattice with 



iV„tp = L,rL„ sites. First there are 2 



2Ns 



choices of 



Sij. We note that there are 2^"*= different Z^ gauge 
transformations. But the constant gauge transformation 
G(%) = — 1 does not change Sjj. Thus there are 2^"'=/2 
different Sjj's in each Zi gauge equivalent class. There- 
fore, there are 2 x 2^"'= different Z^ gauge equivalent 
classes of Sij's. We also note that 

^(_)i.+i.J](_a?Af)(-*A^Af) 



(54) 



different classes of s^j's, 2 



Thus, among the 2 x 2^= 

of them satisfy Hi = (—\^^^+^v and have 

even numbers of -^i^j and V'2,-t fermions. The other 2^='*= 
of them satisfy WiSi^iJ^xSi^i+y = — (— )^^+^b and have 
odd numbers of ipi_i and 'ip2,i fermions. 

For each fixed Sij, there are 2^"*" many-body states 
of the ipi fermions, ie n in \{sij},'i'n) runs from 1 to 
2^site_ Among those 2^"'= many-body states, 2^"'=/2 
of them have even numbers of fermions and 2^"*' /2 
of them have odd numbers of tpi fermions. In order 
for the projection 'P|{sij}, 5'„) to be non-zero, the to- 
tal number of fermions must be even. A physical state 
has even numbers of ("01,11 V'2,i) fermions and even num- 
bers of ^pi fermions, or it has odd numbers of ("01, ij "02,1) 
fermions and odd numbers of ipi fermions. Thus there 
are 2^-'= x 2^"" /2 + 2^"*' x 2^"'= /2 = 4^"*= distinct 
physical states that can be produced by the projection. 
Thus the projection produces all the states in the physi- 
cal Hilbert space. 



B. Physical properties of the spin-ii model 

Let us define a closed-string operator to be 

W{Cclose) - Ui,i,Ui2i,-Ui.il (55) 




FIG. 7; A particle can hop between different sites i,j,k,l. 



where Cdose is an closed oriented string Cdose = ii ^ 
22... in ii formed by nearest neighbor links. Since 
Gciose can intersect with itself, Cdose can also be viewed 
as a closed string-net. We will also call W{Gdose) a 
closed-string-net operator. 

The closed-string-net operators act within the phys- 
ical Hilbert space and commute with the Hamiltonian 
Eq. (41). Thus there is a string- net condensation since 
{W{Gdose)) = ±1 in the ground state of Eq. (41). The 
above strings correspond to the T3 string discussed in 
section IIIC. Unlike the spin-1/2 model, we do not have 
condensed Tl and T2 closed strings in the spin-i-i model. 



We can also define open-string operators 
within the physical Hilbert space 



2 2 
that 



act 



W{Copen) =V4i^ili2'^i2i3---t^i,i-li„V'i„ (56) 

where Gopen is an open oriented string Copen — ii 
22... in formed by nearest neighbor links. W{G) cor- 
respond to the open T3 string defined in section IIIC. 
Just hke the spin-1/2 model Eq. (30), the ends of such 
strings correspond to gapped fermions (if jg] » |t|). The 
ends of W strings only differ from the ends of W strings 
by a local bosonic operator. Thus the ends of W strings 
are also fermions. 

To really prove the ends of W strings are fermions, we 
need to show the hopping of the ends of W strings satisfy 
the fermion hopping algebra introduced in Ref. [14]: 



[tij , tki] —0, if i, j, k, I are all different. 



(57) 



where tji describes the hopping from site i to site j. It 
was shown that the particles are fermions if their hopping 
satisfy the algebra Eq. (57). We note that the ends of 
the W strings live on the sites. The labels in the 

above equation correspond to lattice sites i,j,.... The 
hops between sites i,j,k,l in Fig. 7 are given by 



iipi, a = ±x,±y 



(58) 



Note that the hops between nearest neighbors are taken 
from the Hamiltonian Eq. (41). Since Uij commute with 
each other, the algebra of the above hopping operators 
is just that of fermion hopping operators. In particular, 
the above hopping operators satisfy the fermion hopping 
algebra Eq. (57). Hence, the ends of the W strings are 
fermions. 
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For each fixed configuration Sij, there are 2^"'=/2 
different states (witli even or odd numbers of total ip 
fermions). Their energy are given by the fermion hopping 
Hamiltonian Eq. (43). Let Eo{{sij}) be the ground state 
energy of Eq. (43). The ground state and the ground 
state energy of our spin-ii model Eq. (61) is obtained 
by choosing a configuration Sij that minimize Eo{{sij}). 
We note that Eo{{sij}) is invariant under the Z2 gauge 
transformation Eq. (53). 

When g ^ \t\, the ground state of Eq. (41) has Fi = 
— 1 which minimize the dominating g^^Fi term. The 
ground state configuration is given by 

Si.i+a: = (-)'% Si,i+y = l. (59) 

The il>i fermion hopping Hamiltonian Eq. (43) for the 
above configuration describes fermion hopping with tt- 
flux per plaquette. The fermion spectrum has a form 

Ek = ±2^t^ si-a^ik^) + sm^{ky). (60) 

The low energy excitations of such a hopping Hamilto- 
nian are described by two two-component massless Dirac 
fermions in 2+lD. We see that the ends of the W strings 
are massless Dirac fermions. 

Our model also contain Z2 gauge excitations. The Z2 
vortices are created by flipping Fj = — 1 to Fj = 1 in 
some plaquettes. The Z2 vortex behaves like a 7r-flux to 
the gapless fermions. Thus the gapless fermions carry a 
unit Z2 charge. The low energy effective theory of our 
model is massless Dirac fermions coupled to a Z2 gauge 
field. 



C. Projective symmetry and massless fermions 

We know that symmetry breaking can produce and 
protect gapless Nambu-Goldstone modes. In Ref. [9, 16], 
it was proposed that, in addition to symmetry break- 
ing, quantum order can also produce and protect gapless 
excitations. The gapless excitations produced and pro- 
tected by quantum order can be gapless gauge bosons 
and/or gapless fermions. In this paper we show that 
the quantum orders discussed in Ref. [9, 16] are due to 
string-net condensations. Therefore, more precisely it is 
string-net condensations that produce and protect gap- 
less gauge bosons and/or gapless fermions. The string- 
net condensations and gapless excitations are connected 
in the following way. Let us consider a Hamiltonian that 
has a symmetry described by a symmetry group SG. 
We assume the ground state has a string-net conden- 
sation. Then, the hopping Hamiltonian for the ends of 
condensed string will be invariant under a larger group - 
the projective symmetry group PSG, as discussed in sec- 
tion IV B. PSG is an extension of the symmetry group 
SG, ie PSG contain a normal subgroup IGG such that 
PSG/IGG = SG. The relation between PSG and gap- 
less gauge bosons is simple. Let Q be the maximum 
continuous subgroup of IGG. Then the gapless gauge 
bosons are described by a gauge theory with G as the 
gauge group. [9, 15] Some times the ends of strings are 



fermions. However, the relation between gapless fermions 
and PSG is more complicated. Through a case by case 
study of some PS'G's[9, 16], we find that certain PSG's 
indeed guarantee the existence of gapless fermions. 

In this section, we are going to study a large family 
of exact soluble local bosonic models which depends on 
many continuous parameters. The ground states of the 
local bosonic models have a string-net condensation and 
do not break any symmetry. We will show that the pro- 
jective symmetry of the ends of condensed strings pro- 
tects a massless fermion. As a result, our exact soluble 
model always has massless fermion excitations regard- 
less the value of the continuous parameters (as long as 
they are within a certain range). This puts the results of 
Ref. [9, 16], which were based on mean- field theory, on a 
firmer ground. 

The exact soluble local bosonic models are the spin-^ | 
model 

TT V yx xy yx xy 

i, 

+ E + + ^•^•) (61) 

i 

where 7"^ and 7=*='" are given in Eq. (48) and Eq. (47). 
We will discuss a more general Hamiltonian later. 

The Hamiltonian is not invariant under x —x par- 
ity Px. But it has a a; ^ —x parity symmetry if 
Px is followed by a spin rotation 7^ <-> 7^. That is 
jP^P^Hi^P^P^y = H with 

7P. =7'-^-^ (62) 
Similarly for y —y parity Py, we have 

^p^PyHirtp^Py)-^ = H with 

^^y — ^y 

^Py=^-^ (63) 

In the fermion representation 7p^ and "yp^ generate the 
following transformations 

7p. : K ^ '^f . ^ i^l 

7P. : ^ Af , ^ 4. (64) 

Now let us study how the symmetries T^^y and 
lPmyPx,y are realized in the hopping Hamiltonian 
Eq. (43) for the ends of condensed strings. As discussed 
in section IV B, the hopping Hamiltonian may not be 
invariant under the symmetry transformations Tx^y and 
lPa,,yPx,y directly. The hopping Hamiltonian only has 
a projective symmetry generated by a symmetry trans- 
formation followed by a Z2 gauge transformation G{i). 
Since the 7r-flux configuration does not break any sym- 
metries, we expect the hopping Hamiltonian for the tt- 
flux configuration to be invariant under GxT^, GyTy, 
Gp^lP^Px, and Gp^jp^Py, where Gx,y and Gp^ ^ are the 
corresponding gauge transformations. The action of T^^y 
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and "Yp^^yPx,y on the ijj fermion are given by 



where 





{ix jiy ) 


V'ci^+i.i,,)) 






'^{J'x I'^y ) 


V'ci^.i^+i), 












iPyPy ■■ 






(65) 



For the tt-Aux configuration Eq. (59), we need to choose 
the following Gx,y and Gp^ ^ in order for the combined 
transformation G^yT^^y and Gp^y^p^yP^^y to be the 
symmetries of the hopping Hamiltonian Eq. (43) 



Gx — 1, 

Gp^ = (- 



Gy H-y-' 
Gp, ={-T' 



(66) 



The hopping Hamiltonian is also invariant under a global 
Z2 gauge transformation 



(67) 



The transformations {G(^,Gx,yTx,y,Gp^^^p^ ^Px,y} gen- 
erate the PSG of the hopping Hamiltonian. 

To show that the above PSG protects the masslessness 
of the fermions, we consider a more general Hamiltonian 
by adding 



6Hxi 

2 2 



Y,{t{Ci:i)W{Gii)+h.c) 



(68) 



to Hii 

2 2 



where Cij is an open string connecting site i and 

site j and W{Cij) is given in Eq. (56). The new Hamil- 
tonian is still exactly soluble. We will choose t{Cij) such 
that the new Hamiltonian has the translation symme- 
tries and the Px.y parity symmetries. In the following, 
we would like to show that the new Hamiltonian with 
those symmetries always has massless Dirac fermion ex- 
citations (assuming t(Cy) is not too big comparing to 
9)- 

When t{Cij) is not too large, the ground state is still 
described by the tt-Sux configuration. The new hopping 
Hamiltonian for 7r-flux configuration has a more general 
form 



{ij) 



(69) 



The symmetry of the physical spin-^i Hamiltonian re- 
quires that the above hopping Hamiltonian to be invari- 
ant under the PSG discussed above. Such an invariance 
will guarantee the existence of massless fermions. 
The invariance under GxT^ and GyTy require that 

Xi,i+m = {-y^"'-Xm (70) 

In the momentum space, 

ij 

=eo{k2)5ki-k2 + ^iik2)Ski-k2+Qy (71) 



eo(fe)= 

X —ev 

ei(fc)= ^ e^^ 



rrLx —even 



(72) 



mx=odd 



We note that eo(fe) and ei(fe) are periodic function in the 
Brillouin zone. They also satisfy 

eo(fe) = eo(fe + QJ, ei(fe) = -ei(fe + QJ. (73) 

where = nx and Qy = Try. In the momentum space, 
we can rewrite H as 



(74) 



where — {^k,ipk+Qy)- The sum is over the re- 
duced Brillouin zone: —it < kx < tt and — 7r/2 < ky < 
7r/2. r(fc) has a form 

Note that the transformation7p^ : ^ tp'' changes 

ExijV'lV'j to J2xiji'ii'j with Xij = -~Xji- Thus the 
invariance under Gp^'jp^Px requires that 



-XP.j,P.i = Gp,ii)xijGpx{j) 



or 



X-PxTTi \ ) Am 

In the momentum space, the above becomes 
eo(Pyfc) = - eo(fc) 

ei{Pyk)=ei{k + Qy) (78) 
Similarly, the invariance under Gp^^p^Py requires that 



(76) 



(77) 



-XPyj,Pyi = Gpy{i)XijGpy{j) 



or 



\mxmy -{-niy 



Xr 



X-Pym — -{-Y 

In the momentum space 

eo(^'xfc) 
ei(Pxfc) 

Wc see that the translation Tx.y and x 



eoik + Qy) 
ei(fc) 



(79) 



(80) 



(81) 
'X parity 



iP^Px symmetries of the spin-ii Hamiltonian require 
that eo(fe) = —eo{Pyk) and hence eo(fe)|fej,=o = 0. Sim- 
ilarly, the translation Tx,y and y ^ —y parity iPyPy 
symmetries require that ei(A;)|fe^=o = 0. Thus Tx,y and 
^Px yPx,y symmetries require that r(fe)|fe=o = 0. Using 
Eq.' (73), we find that r(0) = implies that r(Q^) = 0. 
The spin-ii Hamiltonian Eq. (61) has (at least) two two- 
component massless Dirac fermions if it has two trans- 
lation Tx,y and two parity 'yPx,yPx,y symmetries. We 
see that string-net condensation and the associated pro- 
jective symmetry produce and protect massless Dirac 
fermions. 
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VI. MASSLESS FERMIONS AND STRING-NET 
CONDENSATION ON CUBIC LATTICE 

The above calculation and the 2D model can be gen- 
eralized to 3D cubic lattice. We introduce six Majorana 
fermions A", where a = x,x,y,y,z,z. One set of com- 
muting operators on square lattice has a form 



5 \^ \^ 



Using Uij and a complex fcrmion -ipi, we can construct 
an exact soluble interacting Hamiltonian on cubic lattice 



E 

p i a—x,y,z 

^ii ,^2 ,*4 ,*i 



(iUi,i+ai^i'>Pi+a + h.C^j , 



(83) 



where J2p sum over all the square faces of the cubic lat- 
tice, ii, i2, is, and 14 label the four corners of the square 
p. The Hilbert space of system is generated by complex 
fermion operators ijji and 



and there are 16 states per site. 

A physical Hilbert space is defined as a subspacc with 
even numbers of fermions per site. The physical Hilbert 
space has 8 states per site. When restricted in the phys- 
ical Hilbert space, Hi 11 defines our spin-|^| system, 
which is a local bosonic system. 

When g » \t\, our spin-^ii model has two four- 
component masslcss Dirac fermions as its low lying exci- 
tations. The model also has a Z2 gauge excitations and 
the massless Dirac fermions carry unit Z2 gauge charge. 
Again, the model has a string-net condensation in its 
ground state. Both the Z2 gauge excitation and the mass- 
less fermion are produced and protected by the string-net 
condensation and the associated PSG. 



VII. ARTIFICIAL LIGHT AND ARTIFICIAL 
MASSLESS ELECTRON ON CUBIC LATTICE 

In this section, we are going to combine the above 3D 
model and the rotor model discussed in Ref. [44] and 
Rcf. [13] to obtain a quasi-exact soluble local bosonic 
model that contains massless Dirac fermions coupled to 
massless U{1) gauge bosons. 



A. 3D rotor model and artificial light 

A rotor is described by an angular variable 9. The 
angular momentum of 6, denoted as /S^, is quantized as 



integers. The 3D rotor model under consideration has 
one rotor on every link of a cubic lattice. We use ij 
to label the nearest neighbor links, ij and ji label the 

same links. For convenience, wc will define 9ij = ^O^i 



and S, 



-S^-. The 3D rotor Hamiltonian has a form 



Hrotor = E ( E ^li+a. J + 2 E ('^l.i-l-a)^ 

i \ a J i,a 

+ ffi5^cos(^i ) (85) 
p 

(82) Here i = ( 

ix^iy^iz) label the sites of the cubic lattice, 

and a = ±£c,±y, ±z. The ^ sum over all the square 



faces of the cubic lattice, ii, 12, is, and 14 label the four 
corners of the square p. 

When J = 51 = and f7 > 0, the state with all S?j = 
is a ground state. Such a state will be regarded as a state 
with no strings. We can create a string or a string-net 
from the no-string state using the following string (or 
string-net) operator 



(86) 



where C is a string (or a string-net) formed by the nearest 
neighbor links, and Yic i'' ^ product over all the nearest 
neighbor links ij on the string (or string- net). Since the 
closed-string- net operator Wjj(^i^{Cciose) commute with 
(84) Hrotor when J = gi = 0, Wu(i){Cciose) generate a large 



(7(1) 

sc!t of degenerate ground states. The degenerate ground 
states are described by closed string-nets. 

There is another way to generate the degenerate 
groimd states. We note that all the degenerate ground 



states satisfy J2a ^li+a — 0- Let be the common 

eigenstate oi Oij: 0ij\{Oij}) = 0ij\{6ij}). Then the pro- 
jection into J2a ^ii+a = ^ subspacc: "PK^y}) give us a 
degenerate ground state. We note that 



generate a,U{l) gauge transformation 
where 



Oij 



(87) 



(88) 



Thus two U{1) gauge equivalent configurations 6ij and 
6ij give rise to the same projected state 



moij})=rmj}) 



(89) 



We find that the degenerate ground states are described 

by U(l) gauge equivalent classes of 9ij. The degenerate 
ground states also have a U{1) gauge structure. 

When J = but gi ^ 0, the degeneracy in the 
ground states are lifted. One can show that, in this 
case, VliOij}} is an energy eigenstate with an energy 
51 EpCos(^iii2 + 9i^is + 9isi^ + Oi^iJ. Clearly two U{1) 

gauge equivalent configurations 6ij and 9ij have the same 
energy. The non-zero gi makes the closed string-nets to 
fluctuate and vanishing J means that the strings in the 
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string-nets have no tension. Thus J = ground state has 
strong fluctuations of large closed string-nets, and the 
ground state has a closed-string-net condensation. [13] 

When J ^ 0, 'P|{^ij}) is no longer an eigenstate. The 
fluctuations of 6ij describe a dynamical U{1) gauge the- 
ory with 9ij as the gauge potential. [15, 44] 

B. (Quasi-)exact soluble QED on cubic lattice 

To obtain massless Dirac fermions and U{1) gauge 
bosons from a local bosonic model, we mix the spin-^^l 
model and the rotor model to get 

i \ a / i,a 

+ 5i;^cos($p)+5 5^Fp (90) 
p p 

i a=x,y,z 

where $p = Oi^i^ + Oi^i^ + di^i^ + hiii- If we restrict 
ourselves within the physical Hilbcrt space with even 
numbers of fermions per site, the above model is a lo- 
cal bosonic model. 

Let us first set J = 0. In this case, the above model 
can be solved exactly. First let us also set [7 = 0. In this 
case Oij and Uij commute with Hqed and commute with 
each other. Let \{Oij, Sij}, n) be the common eigenstates 
of Oij andX/ij, where n ~ 1, 2, 2^"*= labels different 
degenerate common eigenstates. Within the subspace 
expanded by \{9ij,Sij}, n) , n = 1, 2, 2^^"*-% the Hqed 
reduces to 

Hhop =91 J2 cos(«'p) +9j2fp (^^) 
p p 

i a—x^y,z 

which is a free fermion hopping model. Let 
Sij}, "fn) be the many-fcrmion eigenstate of the 
above fermion hopping model and let E{{6ij,Sij},n) be 
its energy. Then \{Oij), Sij}, ^'„) is also an eigenstate of 
-ffQBD|j=o,i7=o with energy E{{0ij, Sij},n). 
We note that 

Ni = 4i^i + J2Sli+a (92) 

a 

commute with each other and commute with Hqed- 
Thus the eigenstates of Hqed\j=o can be obtained from 

the eigenstates of Hqed\.j=o,u=o by projecting into the 
subspace with Ni = Ni: 

'P{N,}\{Oij,Sij},^n) (93) 

The above state is an eigenstate of Hqed\j=o with an 
energy 

UYNi + E{{eij,Sij},n). (94) 



Eq. (93) and Eq. (94) are our exact solution of 
Hqed\j=o- (We have implicitly assumed that 'P{Ni} also 
perform the projection into the physical Hilbert space of 
even numbers of fermions per site.) 

When U is positive and large, the low energy excita- 
tions only appear in the sector Ni — 0. Those low energy 
eigenstates are given by VHOijjSij}, where P is the 
projection into the Ni = subspace and the even-fermion 
subspace. Their energy is E{{9ij, Sij},n). 

Let us further assume that — gi \t\ and 5 » \t\. In 
this limit, the ground state have fp = —1 and $p = 0. 
We can choose 

6i,i+a =0, a = x,y,z 

Si,i+y =(^) ! 

si,i+z ={-y^+'y (95) 

to describe such a configuration. For such a con- 
figuration, Eq. (91) describes a staggered fermion 
Hamiltonian.[24, 49, 50] The ground state wave func- 
tion PliOij, Sij},'^o) is an eigenstate of the C/(l) closed- 
string-net operator Wu(i){Cciose) with eigenvalue 1. It 
is also a eigenstate of the Z2 closed-string-net operator 
W^(Cciose) with eigenvalue i—)^" where Np is the number 
of the square plaquettes enclosed by Cdose- We see that 
there is a condensation of closed U{1) and Z2 string-nets 
in the J = ground state. In such a string-net condensed 
state, there are gapless fermionic excitations, which are 
described by fermion-hopping in 7r-flux phase. 

In the momentum space, the fermion hopping Hamil- 
tonian Eq. (91) for the 7r-flux configuration has a form 

Hhop =13'*Ur(fc)vl/„,fe + Const. (96) 

k 

where 

Tik) =2tism{k^)Ti + sm{ky)T2 + sin(fc^)r3) 

and Fi = (g) t°, T2 = t'^ 'E) , and F3 = (g) Here 
^1,2,3 g^j.g ^jjg Pauli matrices and t° is the 2 by 2 identity 
matrix. The momentum summation is over a range 
kx e (-7r/2,7r/2), ky G (-7r/2, 7r/2), and k^ € (-7r,7r). 
Since {rjjFj} = 2Sij, i,j = 1,2,3, we find the fermions 
have a dispersion 

E{k) = ±2t^sm^{kx) + sin^{ky) + sm^{k^) (97) 

We see that the dispersion has two nodes at fc = and 
fe = (0,0, tt). Thus, Eq. (91) will give rise to 2 massless 
four-component Dirac fermions in the continuum limit. 

After including the U{1) gauge fluctuations described 
by 9ij and the Z2 gauge fluctuations described by Sij , the 
massless Dirac fermions interact with the U{1) and the 
Z2 gauge fields as fermions with unit charge. Therefore 
the total low energy effective theory of our model is a 
QED with 2 families of Dirac fermions of unit charge 
(plus an extra Z2 gauge filed). We will call those fermions 
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artificial electrons. The continuum effective theory has a 
form 

£ =^^L>o7V/ + Vf^PiDayi 

+ - logiB^ + ... (98) 

where Iq is the lattice constant, / = 1,2, Dq = dt + iao, 
A = di + iai|j=i,2,3, Vf = 2lot, 7''|;:i=o,i,2,3 are 4 X 4 
Dirac matrices, and ■ipi = tjj\'y^. 

We like to pointed out the constant C is Eq. (98) is of 
order 1. Thus the coefficient of the E'^ term oo 

Jlo 

when J = 0. For a finite J, the U{1) gauge field will 
have a non-trivial dynamics. We also like to point out 

that, without fine tuning, the speed of artificial hght, 
Co ~ loVJgij and the speed of artificial electrons, Vf, do 
not have to be the same in our model. Thus the Lorentz 
symmetry is not guaranteed. 

We would like to remark that, for finite J, the U{1) 
closed-string operators no longer condense. A necessary 
(but not sufficient) condition for closed strings to con- 
dense is that the ground state expectation value of the 
closed-string operator satisfy the perimeter law 

(W^[/(i)(Ccio«e)) = Ae-^«/« (99) 

where Lc is the length of the closed string and [A, ^) 
are constants for large closed strings. We note that the 
closed-string operators are the Wilson loop operator of 
the [/(I) gauge field. If the 3-l-lD ?7(1) gauge theory is 
in the Coulomb phase where the artificial light is gapless, 
it was found that[18] 

{WiCiose)) = A(C)e-^^/« (100) 

where A{C) depend on the shape of the closed string 
Cciose even in the large-string limit. Thus the closed 
strings in our model do not exactly condense. The U{1) 
Coiilomb phase is, in some sense, similar to the algebraic- 
long-range order phase of 1-1-1 D interacting boson model 
where the bosons do not exactly condense but the boson 
operator has an algebraic-long-range correlation. 

C. Emerging chiral symmetry from PSG 

Eq. (98) describes the low energy dynamics of the ends 

of open strings (the fermion and the "condensed" 
closed string-nets (the U{1) gauge field). The fermions 
and gauge boson are massless and interact with each 
other. Here we would like to address an important ques- 
tion: after integrating out high energy fermions and 
gauge fluctuations, do fermions and gauge boson remain 
to be massless? In general, interaction between massless 
excitations will generate a mass term for them, unless 
the masslessness is protected by symmetry or some other 
things. We know that due to the U{1) gauge invariance, 
the radiative corrections cannot generate counter term 
that break the U{1) gauge invariance. Thus radiative 
corrections cannot generate mass for the U(l) gauge bo- 
son. For the fermions, if the theory has a chiral symmetry 



tpj e^^^ tpj, 7^ = 7"7"'"7^7'', then the radiative correc- 
tions cannot generate counter terms that break the chiral 
symmetry and thus cannot generate mass for fermions. 
Although the low energy effective theory Eq. (98) ap- 
pears to have the chiral symmetry, in fact it does not. 
This is because that Eq. (98) is derived from a lattice 
model. It contains many other higher order terms sum- 
marized by the ... in Eq. (98). Those higher order terms 
do not have the chiral symmetry. To see this, we note 
that the action of 7^ on '^a,k is realized by a 4 by 4 
matrix 7^ oc rir2r3 oc r^. We also note that the 
periodic boundary conditions of ^'o,fe in the reduced Bril- 
louin zone are given by 

^a,k+Q, = T^®T°*a,fe, *a,fe-KQ„ = T°<E>T^^a,k, (101) 

We find that the action of 7'' is incompatible with the 
periodic boundary conditions since 7^ does not commute 
with (g) r" and r° (g) . Therefore the chiral symmetry 
generated by 7^ cannot be realized on lattice. Due to 
the lack of chiral symmetry, it appears that the radiative 
corrections can generate a mass term 

6£ = iJi,am^i,a (102) 

which is allowed by the symmetry. 

The lack of chiral symmetry on lattice makes it very 
difficult to study massless fermions/quarks in lattice 
gauge theory. In last a few years, this problem was solved 
using the Ginsparg- Wilson relation. [5 1-54] In the follow- 
ing, we would like to show that there is another way 
to solve the massless- fermion/chiral-symmctry problem. 
We will show that our model has an emerging chiral sym- 
metry that appear only at low energies. The low energy 
chiral symmetry comes from the non-trivial quantum or- 
der and the associated PSG in the string-net condensed 
ground state. [9, 15, 16] The Dirac operator in our model 
satisfies a linear relation 

WDW'' =D, We PSG (103) 

in contrast to the non-linear Ginsparg- Wilson relation 

£>7^ + -f^D = aD^^D (104) 

Because of the low energy chiral symmetry, all the 2 fam- 
ilies of Dirac fermions remain massless even after we in- 
clude the radiative corrections from the interaction with 
the U{1) gauge bosons. 

To see how the string-net condensation and the related 
PSG protect the massless fermions, we follow closely the 
discussion in section VC. The Hamiltonian Eq. (90) is a 
mixture of the rotor model and the spin-^^i model. The 
symmetry properties of the rotor part is simple. Here, 
we will concentrate on the spin-i^i part. Eq. (90) is 
not invariant under the six parity transformations: Px.y,z 
and Pxy,yz,zx that generate x ^ —x, y ^ —y, z ^ —z, 
X y, y z, and z x. But it is invariant under par- 
ity Px.y,z and Pxy,yz,zx foUowcd by spin rotations 7p^ ^ 
and 7p^j^ respectively. In the fermion representation 
'yPx,y,z and 7Pj,^_^^,^^ generate the following transforma- 
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tions 



7P« : K 


--Af, 




"fPy ' A- 


-Af, 




IP. ■■ Af 


--A|, 


A ^ V'I, 


7P.. : A- 


-Af, 


Af-Af, 




--A|, 


Af - A|, 


7p : A| 

1 -'^ zx 2. 


-A|, 


A| - Af , 


symmetries 




7Px,«„ 



(105) 

The symmetries T^^^j^, jp^^^^^Px,y,z, and 
7P.y.y.,..Pxv,yz,zx sre realized in the hopping Hamilto- 
nian Eq. (91) through PSG. The hopping Hamihonian 
is invariant only under the symmetry transformations 
followed by proper Z2 gauge transformations G{i). 



Ill 
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sector 



Since the 7r-flux configuration Sij of the spin 
and the zero-flux configuration 6ij of the rotor sector 
do not break any symmetries, we expect the hopping 
Hamiltonian Eq. (91) to be invariant under Gx,y,zTx.y,z, 
Gp.,y_,lP.,y,,Px,y,z, and Gp.^_^ _^^^^jp^y^y^^.^Pxy,yz,zx- 
The action of Tx.y and "fp^^ ^^Pxy,yz,zx on the ip 
fermion are standard coordinate transformations. The 
action of ^p^ y^^Px,y,z on the V fermion are given by 



iP.Px 

iPyPy 
^PyPy 



,iy,iz)' 

-iy,iz)' 



(106) 



For the 7r-flux configuration Eq. (95), we need to choose 
the following Gx,y^z, Gp^^^, and Gp^^^^^^ in order for 
the combined transformation Gx,yTx^y, Gp^ ^jp^ „P^ 
and Gp^^_^^^^7p^^„,^^^P3 



xy,yz,zx 



to be the symmetries of 



the hopping Hamiltonian Eq. (91) 



Gx =(-)'« 
Gp. =(-)^% 



Gy 

GPy 

Gp.. 



<-y 
-i-y 



Gz 
Gpz 
Gp.^ 



-l (107) 

.(^_Yxiy+iviz+izix 



The hopping Hamiltonian is also invariant under a global 
Z2 gauge transformation 



Go : Vi 



-V'i 



(108) 

lP.,y,zPx,y,zy 



The transformations {Gx,y,zTx,y,z, Gp,^ 
Gp,„,„,,„7P.„,„.,„^xy,yz,zx,Go} generate a PSG (a part 
of the full PSG) of the hopping Hamiltonian. 

To study the robustness of massless fermions, we con- 
sider a more general Hamiltonian by adding 



6H = J2 {t{Gij)Wuii){Cij) + h.c) 



(109) 



to Hqed, where Cij is an open string connecting site i 
and site j and W^[/(i)(Gij) an open-string operator 

(110) 

The new Hamiltonian is still exactly soluble, when J = 0. 
We will choose t{Cij) such that the new Hamiltonian has 



the translation symmetries and the Px,y.z parity sym- 
metries. We find that the resulting projective symme- 
try imposes enough constraint on the hopping Hamil- 
tonian for the ends of condensed strings such that the 
Hamiltonian always has massless Dirac fermions (assum- 
ing t{Cij) is not too big comparing to g and gi). Despite 
the PSG transformations Gp^^^^^^^^p^^ ,^^_^^Pxy,yz,zx are 
not needed for the existence of the massless fermions, we 
will still include them in the following discussion. 

For small t{Gij), the ground state is still described by 
the TT-flux configuration. The new hopping Hamiltonian 
for TT-flux configuration has a more general form 



H = Y^iXiM^i + h.c.) 



(111) 



The symmetry of the generalized Hqed requires that 
the above hopping Hamiltonian to be invariant under the 
PSG generated by {Go,Gx,y,zTx.y,z, Gp^ ^^■^p^ ^^Px,y,z, 

Gp^,j^„^.^3.7Px,y.!yz,zi: Pxy,yZ,Zx} ■ 

The invariance under Gx.y.zTx,y,z require that 

Xi,i+m = {-fy--'{-f^^^y+^^)xm (112) 
In the momentum space, 

ij 

= XI ^al3{k2)Ski-k2+aQ^+0Qy (113) 
a,/3=0,l 



where 



eoo(fe) = X] 

my +mz =even^mz =even 

eio(fe) = E 

TJiy +mz =odd,mz =even 

eoi(fe)= X 

my +mz =even^mz =odd 

eii(fe) = X e*-xm. 

'my-\-'m z—odd,m z—odd 



e^'-^^Xm, 



(114) 



We note that eaj3{k) are periodic function in the lattice 
Brillouin zone — tt < kx^y.z < tt. They also satisfy 

e«/3(fe) = {-rec,0{k+Qy+Q,), e„^(fe) = (-)^e„0(fc+Q,, 

(115) 

The r(fe) in Eq. (96) now has a form 



r(fe) = (116) 

/eoo(fe) eio(fe + g^) em{k + Qy) eii(fc + + <5y)\ 

eio(fe) eoo(fe + Q.;) tii{k + Qy) em{k + Qx + Qy) 

eoi{k) eii(fe + g^) eoo{k + Qy) ew{k + Q^ + Qy) 

\eii(fe) eoi(fe + g^) ew{k + Qy) eoo{k + + Qy) J 



Just as discussed in section V C, the invariance under 
Gp^'jP^Px requires that 



-XP.j,P.i = Gp^ii)xijGp,ij) 



(117) 
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or 

X-P.m = +-v-^+™»™«(-)™»Xm (118) 

In the momentum space, the above becomes 

im{-Pxk) = - eoo(fe + Q^)-, 
eio(--Pxfc) = - eio(fe), 
eoi(-Pxfe) =eoi(A; + Q^), 
eii(-P,fe)=-eii(fe). (119) 

where = -kz. Similarly, the invariance under 

GpylPyPy requires that 

X-Pym = -{-)^^^y+^v^^+^'^^{-)^vxm (120) 
In the momentum space 

^ao{- Pyk) = - eoo(fc), 
eio(-Pyfe) =ew{k + Q^), 
eoi{-Pyk) = - eoi{k), 
eni-Pyk)=-eii{k + QJ. 

The invariance under Gp^-yp^Pz requires that 



(121) 



p = ( •^mxmy+mymz+mzmx ^ )^'X (1^2) 

In the momentum space 

€oo{-Pzk) = - eoo(fc), 
eio(-^'^fe) =-eio(fc + Q^), 
eoi(-Pjfe) = - eoi(A;), 
eii(-P.fe)=eii(fe + QJ. (123) 

The invariance under Gp^^jp^^Pxy requires that 

XPxyi,Pxy3 = Gp^y {i)XijGp^y (j) (124) 



or 



XPtcym = {-T 

In the momentum space 



(125) 



eoo{Pxyk) =eoo(fe), 
eio{Pxyk) =eio(fe + QJ, 

eOl{Pxyk) =£01 (fe + QJ, 

£n{Pxyk) =eii{k). 



(126) 



The invariance under Gp^^jp^^Pyz requires that 

XPy.m = i-r^'^'Xm (127) 

or 



eoo{Pyzk) =eoo{k), 
eio{Pyzk) = - eio(fc), 
eoi{Pyzk) = - eoi(fe), 
eii{Pyzk) =eii(fe). 



The invariance under Gp^^^p^^Pzx requires that 

XP.xm = (-)™»™.+™v™^+™^'"^Xm (129) 



or 



eoo{Pzxk) =£00 (fe), 
eioiPzxk) =£io(fc + QJ, 
ioi{Pzxk) = - £oi(A;), 
euiPzxk) =cn{k + Q^). 



(130) 



We see that Eq. (119) require that €io{k)\ky=k^=o = 
and eii{k)\ky=k^=o = 0. Eq. (123) require that 
eoo(fe)|fc,=fe„=o = and £oi(fe)|fc^=fc^=o = 0. Thus the 
e„^(0) = 0. When combined with Eq. (115), Eq. (119), 
and Eq. (123), we find 

eaf3{axQx+ayQy+azQz) =0, a^, ay, az = 0,1 (131) 

Therefore r(fe) = when fe = 0, Q^,. The two transla- 
tion Tx^y and the three parity lPx,y,zPx,y,z symmetries 
of Hqed guarantee the existence of at least 2 four- 
component masslcss Dirac fcrmions. Or more precisely, 
no symmetric local perturbations in the local bosonic 
model Hqed can generate mass terms for the two mass- 
less Dirac fcrmions in the impcrturbcd Hamiltonian. 

Since the mass term in the continuum effective field 
theory is not allowed by the underlying lattice PSG, we 
say that our model has an emerging chiral symmetry. 
The masslessness of the Dirac fermion is protected by 
the quantum order and the associated PSG. 



VIII. QED AND QCD FROM A BOSONIC 
MODEL ON CUBIC LATTICE 



In this section, we arc going to generalize the results 
in Ref. [55] and Ref. [15] and use a bosonic model on 
cubic lattice to generate QED and QCD with 2Nf fam- 
ilies of massless quarks and leptons. To describe the 
local Hilbcrt space on site i in our bosonic model, it 
is convenient to introduce fermions A? and V"?"' where 
o = 1, Nf, n = l, 2Nf and a = 1, 2, 3. A? is in the 
fundamental representation of a SU{Nf) group. V-""" is 
in the fundamental representation of a SU (3) color group 
and a SU{2Nf) group. The Hilbert space of fermions is 
bigger than the Hilbert space of our boson model. Only a 
physical subspace of the fermions Hilbert space becomes 
the Hilbert space of our boson model. The physical states 
on each site is formed by color singlet states that satisfy 



AfA?j"^' + vrv: 



S"'^Nf]\%ny,)=0 (132) 



(128) 



where Nj is assumed to be even. Once restricted within 
the physical Hilbert space, the fermion model becomes 
our local bosonic model. 

In the fermion representation, the local physical oper- 
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ators in our bosonic model are given by 



amn _^jjna\ i na _ _}_xmn ila] jla 



(133) 



We note that by definition M?° = 5?" = 0. The Hamil- 
tonian of our boson model is given by 

^ =^ E ^r^r + ^ E '^rf 



(ij> 



iii) 



Let us assume, for the time being, J3 = 0. In terms of 
fermions, the above Hamiltonian can be rewritten as 

^ = - # E r/r,'''i^r^7'' - # E m'>^''i>^1 

^ {ij) ^ (id) 

+ Const. (135) 
Using path integral, we can rewrite the above model as 

Z=J I?(V/)2?(V^)I?(ao)I?(«)2?(x)e^^*^ 

(iJ) 

{ii) 

-f E^KH)-f Exi.-Xl- (136) 

(ij) (ij) 

where (V-"?)^ — i^T^ ^'^T'^ ^'^T^)^ ^'^^ flo(*) and are 
3x3 complex matrices that satisfy 



Aj = '^ji^ «o(i) = al{i) 



(137) 



When J3 ^ 0, the Lagrangian may contain terms that 
mix Xij and Uij: 

L =i;fi[dt + iaoUMf - Y^i^fuijij] + h.c.) 

(ij) 

+ \fi[dt + i'n:ao(i)]A? - Y.^\f^ij^'j + ^^•^•) 

(ij) 

-^Y.^(uii4j)-^Y.xijx\i 

+ ^[Xi,- det(w,-i) + h.c.\ (138) 

where C is a 0(1) constant. We note that the above 
Lagrangian describes a U{1) x SU{3) lattice gauge theory 
coupled to fermions. 

The field ao{i) in the Lagrangian is introduced to en- 
force the constraint 

^na^nff _ ^na\ ^ _ A?Af ^"'^ = (139) 



As in standard gauge theory, the above constraint really 
means a constraint on physical states, ie all physical 
states must satisfy 



XfKd'^P + VrVf - 5'^\Ns ) 1$^,,,) = (140) 



The above is the needed constraint to obtain the Hilbert 

space of our bosonic model. 

Here we would like to stress that writing a bosonic 
model in terms of gauge theory does not imply the ex- 
istence physical gauge bosons at low energy. Using pro- 
jective construction, we can write any model in terms 
of a gauge theory of any gauge group. [28, 56] The exis- 
tence of low energy gauge fluctuations is a property of 
ground state. It has nothing to do with how we write the 
Hamiltonian in terms of this or that gauge theory. 

Certainly, if the ground state is known to have certain 
gauge fiuctuations, then writing Hamiltonian in term of 
a particular gauge theory that happen to have the same 
gauge group will help us to derive the low energy effective 
theory. Even when we do not know the low energy gauge 
fluctuations in the ground state, wc can still try to write 
the Hamiltonian in a form that contains certain gauge 
theory and try to derive the low energy effective gauge 
theory. Most of the times, we find the gauge fluctuations 
in the low energy effective theory is so strong that the 
gauge theory is in the confining phase. This indicates 
that we have chosen a wrong form of Hamiltonian. How- 
ever, if we are lucky to choose the right form of Hamilto- 
nian with right gauge group, then the gauge fluctuations 
in the low energy effective theory will be weak and the 
gauge fields ao, Xij and Uij will be almost like classical 
fields. In this case, wc can say that the ground state of 
the Hamiltonian contains low energy gauge fluctuations 
described by ao, Xij and Uij. In the following, we will 
show that the U{1) x SU{^) fermion model Eq. (136) is 
the right form for us to write the Hamiltonian Eq. (134) 
of our bosonic model. 

After integrating out the fermions, we obtain the fol- 
lowing effective theory for ao{i), Xij and Uij 



V{ao)V{u)e 



J dtNjL^fj{u,ao) 



(141) 



where Lejj does not depend on Nf. We see that, in the 
large TV/ limit, Xiji '^ij and oq indeed becomes classical 
fields with weak fiuctuations. 

In the semi-classical limit, the ground state of the sys- 
tem is given by the ansatz {Xij^Uij,ao{i)) that minimize 
the energy —Lefj. We will assume that such an ansatz 
have TT flux on every plaquette and takes a form 



/(i.i+x ^X: 

Xi.i+s = - «(-)'"+'''x, 



u. 



lU, 



Xi,i+y 



^i,i-\-y 



u. 



z,z-\-z 



i(-y^+^yu, ao{i) =0. 



(142) 



(If the TT-flux ansatz does not minimize the energy, we 
can always modify the Hamiltonian of our bosonic model 
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to make the tt-Aux ansatz to have the minimal energy.) 
Despite the i dependence, the above ansatz actually de- 
scribe translation, rotation, parity, and charge conjuga- 
tion symmetric states. This is because the symmetry 
transformed ansatz, although not equal to the original 
ansatz, is gauge equivalent to the original ansatz. 
The mean-field Hamiltonian for fermions has a form 

H = Y^i^pfuij^PJ + XfxijX^j + h.c.) (143) 

{ij) 

The fermion dispersion has two nodes at A; = and 
k = (0,0,7r). Thus there are 2Nf x 7 massless four- 
component Dirac fermions in the continuum limit. They 
correspond to quarks and leptons of 2Nf different fami- 
lies. Each family contains six quarks (two flavors times 
three colors) that carry SU{3) colors and charge 1/3 for 
the U{1) gauge field, and one lepton that carry no SU{3) 
colors and charge 1 for the U{1) gauge field. 

Including the collective fluctuations of the ansatz, the 
f/(l) X SU{3) = U{3) fermion theory has a form 

L = i(at + iao{i)WJ + J2 ^fuae^^^^i^] 

i ij 

+ E + ^Trao(i))A,^ + ^ Af xy det(e--)A5 

i ij 

(144) 

where aij are 3x3 hermitian matrices, describing U{1) 
and SU{3) gauge fields. In the continuum limit, the 

above becomes 

+ XlMD'„J°Xl,a + v'f\l,aD[fXl,a (145) 

with vf ~ ^0-^1, v'f ~ hJ2-, -D/i = 9^ + ia^l■, D'^j^ = + 
iTra^, I = 1,2, and are 4 x 4 Dirac matrices. [24, 49, 
50] Xj a and ?/'/.« arc Dirac fermion fields, tpi^n form a 
fundamental representation of color SU{3). 

If we integrate out ao and aij in Eq. (144) first, we will 
recover the bosonic model Eq. (134). If we integrating 
out the high energy fermions first, the U{1) x SU (3) gauge 
field will acquire a dynamics. We obtain the following 
low energy effective theory in continuum limit 

+ Xl.aD'ol^Xj^a + v'jXj^aD'ij'-Xl^a 

+ -^[TrFoiF°' + clTvFijF'^ + ... (146) 

where the velocity of the U{3) gauge bosons is Ca ~ hJi,2, 
and ... represents higher derivative terms and the cou- 
pling constant as is of order 1/Nf. 

In the large N f limit, fiuctuations of the gauge fields 
are weak. The model Eq. (146) describes a [/(I) x SU{3) 
gauge theory coupled weakly to 2Nf families of mass- 
less fermions. Therefore, our bosonic model can generate 
massless artificial quarks and artificial leptons that cou- 
ple to artificial light and artificial gluons. As discussed 
in Ref. [15], the PSG of the ansatz Eq. (142) protects 
the masslessness of the artificial quarks and the artificial 
leptons. Our model has an emerging chiral symmetry. 



IX. CONCLUSION 

In this paper, we studied a new class of ordered states 
- string-net condensed states - in local bosonic models. 
The new kind of orders does not break any symmetry 
and cannot be described by Landau's symmetry break- 
ing theory. We show that different string-net condensa- 
tion can be characterized (and hopefully, classified) by 
the projective symmetry in the hopping Hamiltonian for 
ends of condensed strings. Similar to symmetry breaking 
states (or "particle" condensed states), string-net con- 
densed states can also produce and protect gapless exci- 
tations. However, unlike symmc^try breaking states which 
can only produce and protect gapless scaler bosons (or 
Nambu-Goldstone modes), string- net condensed states 
can produce and protect gapless gauge bosons and gap- 
less fermions. It is amazing to see that gapless fermions 
can even appear in local bosonic models. 

Motivated by the above results, we propose the fol- 
lowing locality principle: The fundamental theory for 
our universe is a local bosonic model. Using several lo- 
cal bosonic models as examples, we try to argue that 
the locality principle is not obviously wrong, if we as- 
sume that there is a string-net condensation in our vac- 
uum. The string-net condensation can naturally pro- 
duce and protect massless photons (as well as gluons) 
and (nearly) massless electrons/quarks. However, to re- 
ally prove the string-net condensation in our vacuum, we 
need to show that string-net condensation can generate 
chiral fermions. Also, the above locality principle has not 
taken quantum gravity into account. It may need to be 
generalized to include quantum gravity. In any case, we 
can say that we have a plausible understanding where 
light and fermions come from. The existence of light 
and fermions is no longer mysterious once we realize that 
they can come from local bosonic models via string-net 
condensations. 

The string-net condensation and the associated PSG 
also provide a new solution to the chiral symmetry and 
the fermion mass problem in lattice QED and lattice 
QCD. We show that the symmetry of a lattice bosonic 
model leads to PSG of the hopping Hamiltonian for the 
ends of condensed strings. If the ends of condensed 
strings are fermions, then PSG can some times protect 
the masslessness of the fermions, even though the chiral 
symmetry in the continuum limit cannot be generalized 
to the lattice. Thus PSG can lead to an emerging chiral 
symmetry that protect massless Dirac fermions. 

In this paper, we have been stressing that string-net 
condensation and the associated PSG can protect the 
masslessness of fermions. However, most fermions in na- 
ture do have masses, although very small comparing to 
the Planck mass. One may wonder where those small 
masses come from. Here we would like to point out that 
the PSG argument for masslessness only works for ra- 
diative corrections. In other words, the fermions pro- 
tected by string-net condensation and PSG cannot gain 
any mass from additive radiative corrections caused by 
high energy fluctuations. However, if the model has in- 
frared divergence, then infrared divergence can give the 
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wo\ild-bc-massless fermions some masses. The acqiiired 
masses should have the scale of the infrared divergence. 
The 3+lD QED model studied in this paper do not have 
any infrared divergence. Thus, the artificial electrons in 
the model are exactly massless. But in the bosonic model 
discussed in section VIII, the SU{3) gauge coupling as 
runs as 

^^5' ^ 11 - (2/3)(2jV/) 
d\n{M^) 477 ^ ' 

where M is the cut-off scale. Thus when Nf < 8, as 
has a logarithmic infrared divergence. In general, for 

models with U{1) and SU{3) gauge interactions and 
a right content of fermions, the SU{3) gauge interac- 
tions can have a weak logarithmic infrared divergence in 
3+lD.[57, 58] This weak divergence could generate mass 
of order e^'^'^^^-'^^^Mp, where Mp is the Planck mass 
or the GUT scale (the cut-off scale of the lattice theory), 
C = 0(1) and as{Mp) is the dimensionless gauge cou- 
pling constant at the Planck scale. An C/as{Mp) ~ 40 
can produce a desired separation between the Planck 



mass/GUT scale and the masses of observed fermions. It 
is interesting to see that, in order to use string-net con- 
densation picture to explain the origin of gauge bosons 
and nearly massless fermions, it is important to have a 
four dimensional space-time. When space-time has five 
or more dimensions, the gauge-fermion interaction do not 
have any infrared divergence. In this case, if a string-net 
condensation produces massless fermions, those fermions 
will remain to be massless down to zero energy. In 2-l-lD, 
the gauge interaction between massless fermions is so 
strong that one cannot have fermionic quasiparticles at 
low energies. [59-61] It is amazing to see that 3-1-1 is the 
only space-time dimension that the gauge bosons and 
fermions produced by string-net condensation have weak 
enough interaction so that they can be identified at low 
energies and, at same time, have strong enough interac- 
tion to have a rich non-trivial structure at low energies. 
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